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blockaiio  and  skin-frinlion  roduction  due  lo  mass  injcclion. 

I he  classie  Crocco  lempcralun— velociU  relation  was  examined  at  liyper- 
sonic  flow  conditions  on  a sharp  porous  com*.  IVepartures  from  the  linear 
Crooco  relation  were  observed  tluit  tended  to  follow  a nuHlified  qnailralic  distri- 
bution in  velocit\ . In  ticneral,  the  data  tended  to  follow  the  linear  C'rocco  re- 
lation willi  a tempcr.'iture  overslioot  resulting  from  the  kinetic  energ\  of  the  flow 
whicli  is  typical  of  flat-plate  type  flows  previoush  investigated,  A non-similarity 
parameter  was  introduceil  to  account  for  the  effect.-,  of  mass  transfer,  pressure 
gradient  and  temperature  p»itential. 

An  evaluation  of  the  boundar.\  la\e.r  profile  shape  (in  terms  of  the  velocity 
power-law  e.\ponent  a-..l  \iseou.s  laser  integral  thickness  relations)  resulted  in 
correlations  for  viscous  la\er  thicknesses  w'hich  account  for  the  effects  of  Mach 
number,  wall  tomperature.  mass  transfer,  injectant  molecular  weight  and  angle 
of  attack. 

flu;  effect  of  Mach  manber.  wall  temperature  and  injtxjtant  molecular  weight 
on  heat  l)lockage  and  skin  friction  re<iuetion  rc.sulting  from  mass  injection  was 
evaluated,  resulting  in  seini-emi)irical  correlations  which  are  based  on  a modifi- 
cation of  a coinpressibilitv  transformation  tlu!or\  expanded  to  incluile  the  latest 
hypersonic  data,  file  re.sultant  heat  blockage  and  skin  fiiction  reduction  corre- 
lations were  verified  b\  comparison  with  data  o\er  a wide  range  of  .Mach  number, 
wall  temixjraturo  and  injectant  molecular  weight. 

The  correlation  dcvclo|)ment.  ba.sed  on  an  analytically  derived  compressi- 
bilit\  transformation  and  an  empirically  derived  vi.scous  transformation,  em- 
plo\e<l  a critical  blowawat  parameter  concept  to  provide  a more  tractable  set 
of  correlation.s  for  engineering  design.  .Angle  of  attack  heat  tr.insfer  was  satis- 
factorily jiredu  ted  using  an  e.sisting  . emi-empirical  turluilcit  l)oundar\  laver 
code  which  utilized  Ihe  aloremenlioned  mass  ir.insfer  coi’relations. 
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FOREWORD 


This  report  presents  a detailed  examination  of  engineering  design  relationships 
for  turbulent  boundary  layer  flows  which  account  for  the  effects  of  mass  transfer,  angle 
of  attack,  and  Reynolds  number.  In  particular,  the  study  is  concerned  with  turbulent 
boundary  layer  characteristics  associated  with  energy  profile  shape  (Crocco  tempera- 
ture-velocity relation),  velocity  power  law  exponent,  viscous  layer  thickness,  skin 
friction  reduction,  and  heat  transfer  blockage.  These  tasks  were  selected  as  a conse- 
quence of  their  direct  use  in  engineering  design  codes.  The  data  base  for  this  study 
consists  of  the  experimental  work  obtained  by  GE-RESD  under  SAMSO  Air  Force  con- 
tract F04701-70-C-0179  entitled  "Strategic  Re-entry  Technology  Program  - Phase  II" 
(STREET-G)  Task  4, 11  (1972)  and  Task  4. 14  (1973)  as  well  as  data  from  the  literature. 
The  work  supported  herein  was  performed  at  the  General  Electric  Company,  Re-entry 
and  Environmental  Systems  Division,  Valley  Forge,  Pennsylvania. 

This  document  represents  the  culmination  of  work  performed  under  Contract  Num- 
ber F33615-74-C-3016,  Project  1366;  '"Aerodynamic  Empirical  Techniques  for  Strategic 
Missiles. " The  work  was  performed  in  the  time  period  November  1973  through  April 
1975.  This  investigation  was  conducted  for  the  United  States  Air  Force  Flight  Dynam- 
ics Laboratory  under  the  technical  direction  of  Mr.  Richard  R.  Smith  (AFFDL/FXG). 

The  final  report  was  submitted  by  the  authors  for  publication  approved  on  April  25,  1975. 

The  authors  wish  tj  acknowledge  the  following  General  Electric  personnel  for  their 
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David  Nestler  and  Ms.  C,  Dougherty  for  computational  assistance  and  finally  to  Mr. 

R.  Sweeney  and  Mr.  W.  Norman  for  computation  and  plotting. 
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TECHNICAL  SUMMARY 

The  pr<Alems  examined  in  this  study  are  concerned  with  turbulent  boundary  layer 
chan«teristics  that  are  required  for  engineering  design  codes.  These  consist  of: 

(1)  bi>undary  layer  profile  shape  (velocity  power  law  esqnxient),  (2)  energy  profile 
shape  (Crocco  temperature  - velocity  relation),  (3)  viscous  layer  thicknesses  (6,6’*, 
and  0 ),  (4)  skin  friction  reducti<m,  ai^  (5)  heat  transfer  blockage.  The  above  were 
investigated  relative  to  Improving  existing  engineering  design  capability  for  turbulent 
boundary  layer  flows  that  are  associated  with  the  effects  of  angle  of  attack,  Re3molds 
number,  compressibility,  and  mass  transfer.  This  study  serves  to  provide  two  prim- 
ary ftmctions:  first,  to  upgrade  industrial  design  codes  that  rely  upon  two-dimensional, 
low  Mach  niunber  data  cou'elations  and  secondly,  to  provide  an  insight  to  sophisticated 
3-D  finite  difference  or  integral  viscous  boundary  layer  codes. 

The  classic  Crocco  temperature-velocity  relation  was  examined  at  hypersonic  flow 
conditions  on  a sharp,  porous  cone.  The  evaluation  was  made  with  the  effects  of  mass 
transfer,  angle -of-attack,  and  molecular  weight  variations  of  the  injectant  gas  for  a 
free  stream  Mach  number  of  approximately  8. 0.  Departures  from  the  linear  Crocco 
relation  were  observed  that  tended  to  follow  a modified  quadratic  distribution  in  velocity. 
In  general,  the  data  tended  tu  follow  the  linear  Crocco  relation  with  a temperature  over- 
shoot resulting  from  the  kinetic  energy  of  the  flow  which  is  typical  of  flat  plate  type  flows 
previously  investigated. 

Specifically,  no  apparent  effect  of  Mach  number  or  Reynolds  number  was  noted. 
Mass  transfer  tended  to  depart  from  the  linear  relationship  toward  a modified  quadratic 
velocity  relation.  The  effect  of  boundary  layer  transition  (up-stream  history)  should  be 
considered  when  evaluating  the  Crocco  variables.  For  windward  angle  of  attack,  the 
data  tended  to  follow  the  linear  distribution;  whereas,  for  leeward  conditions,  the  data 
were  characterized  by  a quadratic  behavior.  Moreover,  for  laminar  boundary  layers 
the  data  was  linearly  oriented,  and  for  turbulent  flow  conditions  the  data  appeared  to 
follow  the  quadratic  relation.  No  specific  conclusions  could  be  made  regarding  the 
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Crocco  relation  for  non-air  injection  cases.  Finally,  a ncm-similarity  parameter  is 
introduced  to  account  for  the  effects  of  blowing,  pressure  gradient,  and  temperature 
potentiaL 

This  portion  of  the  study  is  considered  an  important  contribution  to  boundary  layer 
behavior  inasmuch  as  departures  from  similarity  of  the  thermal  and  hydrod^mamlc 
(Pr  ^ unity)  characteristics  of  the  boundary  layer  can  be  accounted  for  with  slight 
modifications  to  the  linear  Crocco  relation.  In  particular,  transformatlcn  ftmctlons 
which  relate  incompressible  to  compressible  theory  can  be  modified  to  account  for  de- 
partures from  a unit  Prandtl  number  that  include  longitudinal  pressure  gradient,  mass 
turmsfer  and  temperature  potential. 

The  turbulent  boundary  layer  characteristic  parameters  6,  5 $,  and  n were 

also  examined  relative  to  the  influence  of  Reynolds  number,  angle  of  attack,  mass 
transfer,  and  molecular  weight  variations  of  the  injectant  gases.  The  study  empha- 
sized recently  obtained  data  at  hypersonic  flow  conditions  on  a sharp  porous  cone.  In 
general,  it  was  found  that  the  state  of  boundary  layer  development  was  sensitive  to  the 
velocity  power  law  exponent  as  well  as  the  viscous  layer  thicknesses. 

Specifically,  the  power  law  exponent  showed  significant  variations  from  the  classic 

ftilly-developed turbulent  flow  value  of  n = 7.  Asa  consequence,  reference  momentum 

thickness  relations  (i.  e. , n = x = 0)  had  to  be  modified  to  include  a generalized  velocity 

power  law  exponent.  Prediction  capability  for  this  reference  situation  is  shown  to  be 

in  excellent  agreement  with  data.  Concerning  mass  transfer  effects,  correlations  were 

developed  for  the  viscous  layer  thicknesses  which  include  molecular  weight  effects. 

When  combining  angle  of  attack  effects  with  blowing,  a hyperbolic  type  correlation  was 

developed  that  considers  the  viscous  layer  thickness  (when  normalized  by  its  a = 0 

value)  as  a function  of  a/e  as  well  as  the  blowing  parameter. 

c 

The  development  of  the  viscous  layer  thickness  correlations  is  considered  a sig- 
nificant advancement  in  the  state  of  the  art  for  turbulent  flow  engineering  prediction 
capability.  In  particular,  the  ability  to  determine  viscous  layer  parameters  for  three- 
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dimensional  configurations,  includiiig  the  effect  of  mass  transfer  and  molecular  weight, 
affords  the  design  engineer  greater  latitude  than  previously  attained.  Moreover,  these 
correlations  also  allow  for  estimates  of  turbulent  flow  behavior  over  complex  geomet- 
ries such  as  control  surfaces  and  surface  disc<mtinuities. 

Finally,  the  effect  of  Mach  number,  wall  temperature  and  injectant  molecular 
weight  on  heat  blockage  and  skin  friction  reduction  resulting  from  mass  injection  was 
evaluated  for  turbulent  flow  over  surfaces  with  negligible  axial  pressure  gradients. 

The  evaluation  resulted  in  semi-empirical  correlations  which  are  based  on  a modifica- 
tion of  a compressibility  transformation  theory  expanded  to  include  the  latest  hypersonic 
data  up  to  a free  stream  Mach  number  of  8. 0.  The  resultant  heat  blockage  and  skin 
friction  reduction  correlations  Wei's  verified  by  comparison  with  data  over  a wide  range 
of  Mach  number,  wall  temperature  and  injectant  tiolecular  weif^t.  Emphasis  was 
placed  on  recently  acquired  heat  transfer  data  at  hypersonic  flow  conditions  on  a porous 
sharp  cone. 

The  correlation  development,  based  on  an  analytically  derived  compressibility 
transformation  and  an  empirically  derived  viscous  transformation,  employed  a critical 
blowaway  parameter  concept  to  provide  a more  tractable  set  of  correlations  for  engi- 
neering design.  Angle  of  attack  heat  transfer  was  satisfactorily  predicted  using  an 
existing  semi-empirical  turbulent  boundary  layer  code  which  utilized  the  aforemen- 
tioned mass  transfer  correiations. 

By  accounting  for  the  effects  of  Mach  number,  wall  temperature  and  injectant  mol- 
ecular weight  on  heat  transfer  and  skin  friction  attenuation  due  to  mass  transfer,  these 
correlations  contribute  significantly  to  the  state  of  the  art,  effectively  providing  im- 
proved engineering  design  relationships  suitable  fur  inclusion  in  existing  semi-empirical 
turbulent  boundary  layer  codes.  In  addition,  the  theoretical  compressibility  transfor- 
mation concept  employed  in  this  development  has  been  shown  to  be  effective  and  provides 
a sound  fundamental  base  for  future  turbulent  boundary  layer  studies. 
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SECTION  1 
INTRODUCTION 

The  alleviation  of  the  high  heating  rates  encountered  hy  surfaces  of  hypersonic 
vehicles  has  been  recognized  as  an  important  problem.  One  of  the  cooling  methods 
that  has  shown  promise  is  mass-transfer  cooling,  wherein  a "foreign"  material  is 
transferred  from  the  vehicle  surface  into  the  boundary  layer.  This  has  a two-fold 
advantage  in  attenuating  the  heat-transfer  problem.  The  transferred  coolant  may  ab- 
sorb heat  from  the  boundary  layer  through  a phase  change  (sublimation,  evaporation, 
melting,  etc. ) as  well  as  providing  high  thermal  heat  capacitance.  In  addition,  it  has 
been  shown  that  the  introduction  of  a material  (with  its  normal  velocity  component)  at 
the  surface  acts  to  decelerate  the  flow  and,  consequently^  to  reduce  the  skin  friction 
which,  in  turn,  reduces  the  heat  transfer  at  the  walL  The  problem  of  coupling  the 
fluid  mechanics  of  boundary  layer  flow  with  material  response  is  further  enhanced  when 
considering  angle  of  attack  effects.  The  evaluation  of  mass  addition  on  the  aerodynamic 
characteristics  of  re-entry  vehicles  at  angle  of  attack  is  a complicated  but  definable 
process.  Since  the  evaluation  of  viscous  effects  utilizing  classical  boundary  layer  the- 
ories requires  knowledge  of  the  inviscid  flow  conditions  at  the  outer  edge  of  the  bound- 
ary layer,  the  inviscid  streamline  pattern  of  the  flow  and  the  resulting  outer  edge  flow 
conditions  must  be  determined.  The  angle  of  attack  of  the  vehicle  introduces  local  pres- 
sure gradients  and  cross  flow  effects  which  must  be  factored  into  the  analysis.  It  is 
these  effects  which  tend  to  dominate  the  flow  computational  procedures. 

Within  the  aerospace  community,  very  sophisticated  computer  programs  have  been 
developed  (or  detailed  computation  of  viscous  flow  properties.  These  programs  pro- 
vide rigorous  solutions  to  many  of  the  complex  phenomena  encountered  by  a re-entering 
vehicle,  e.  g.,  blunt  body  flows,  separated  flowL  ablation,  etc.  However,  they  usually 
require  significant  computational  time  and  can  be  sensitive  to  input  values  specified  for 
configuration  and  flow  conditions.  Consequently,  other  classes  of  analytical  methods 
have  been  developed.  These  techniques  allow  more  rapid  calculations  for  use  in  design 
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tradeoff  studies,  trajectory  calculations,  or  evaluation  of  experimental  data.  These 
rapid  design  techniques  (engineering  type)  make  use  of  the  results  of  the  more  sophis- 
ticated analyses  and  of  experimental  results  either  directly  or  in  the  form  of  correla- 
tions, to  provide  approximate  but  sufficiently  accurate  results.  It  is  the  engineering 
type  of  solution  that  will  be  examined  in  this  study. 


For  turbulent  flows,  these  engineering  codes  consider  the  heat  transfer  blockage, 
skin  friction  reduction,  and  boundary  layer  thickness  growth  with  mass  addition  and 
angle  of  attack  by  semi-empirical  means.  In  addition,  it  is  customarily  assumed  that 
the  velocity  profile  exponent  (n)  is  7 and  that  the  Crocco  integral  is  valid.  Concern- 
ing the  former,  data  obtained  in  the  SAMSO  sponsored  Strategic  Re-entry  Technology 
(STREET-G)  Program  (corroborated  in  several  other  reports)  indicate  that  a nonsimi- 
lar boundary  layer  exists  downstream  of  transition  (which  could  persist  for  50  to  100 
boundary  layer  thicknesses).  In  this  region,  which  always  exists  on  flight  vehicles, 
there  is  a paucity  of  data.  Furthermore,  no  data  containing  detailed  boundary  layer 
Information  for  ablating  bodies  exists  for  the  angle  of  attack  case.  On  the  other  hand, 
the  classic  Crocco  temperature -velocily  relationship  has  been  traditionally  used  to 
relate  the  density  variations  in  a boundary  layer  to  the  velocity.  In  so  doing,  one  must 
postulate  an  equation  of  state  together  with  the  boundary  layer  assumption  for  zero 
pressure  gradient  normal  to  the  flow  direction.  One  such  method  of  relating  the  den- 
sity and  velocity  is  to  assume  unit  Prandtl  and  Schmidt  numbers  throughout  the  bound- 
ary layer.  Often  referred  to  as  the  Reimolds  Analogy,  this  assumption  implies  a re- 
lationship exists  between  shear  stress,  heat  transfer,  and  mass  transfer.  Moreover, 
the  Crocco  relationship  was  developed  from  the  two-dimensional  boundary  layer  equa- 
tions and  its  use  in  three-dimensional  flows  and  flows  with  mass  transfer  must  be  care- 
fulh  examined  to  determine  its  range  of  applicability  for  engineering  prediction  tech- 
niques. 

As  indicated  above,  the  emphasis  of  this  study  is  to  upgrade  and  supplement  data  cor- 
relations that  are  required  in  engineering  design  codes.  In  particular,  the  recent  hyper- 
sonic mean  flow  and  surface  data  reported  in  References  1 and  2 will  be  emphasized. 
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The  STREET -G  experiment  was  performed  at  the  Arnold  Engineering  Development 
Center  in  Tunnel  B at  Mach  8 using  a five-degree  half  angle  porous  cone  which  is  five 
feet  long.  The  experiment  considered  mass  addition  rates  which  cover  the  ablation 
gas  range  for  current  heat  shield  materials  (simulated  ablation  via  gaseous  mass  injec- 
tion of  nitrogen,  argon,  Freon  12,  and  helium  through  the  porous  surface).  Other 
significant  parameters  that  were  varied  include  molecular  weight  of  the  injectant, 
Reynolds  number,  and  vehicle  angle  of  attack  (Of  to  lOf ).  These  data  are  considered 
unique  and  represent  an  advancement  m the  state  of  the  art  for  turbulent  boundary  layer 
characteristics. 

The  objectives  of  the  proposed  study  are  to  provide  improved  engineering  design 
relationships  for  turbulent  boundary  layer  flows  to  account  for  the  effects  of  mass  trans- 
fer, angle  of  attack,  and  Reynolds  number.  Several  correlations  employing  the  effects 
of  ablation  (mass  transfer)  on  local  viscous  flow  properties  have  been  reported  in  the 
literature.  These  correlations  are,  in  general,  based  on  incompressible  "flat  plate" 
type  theory  together  with  low  and  high  speed  data  (supersonic).  In  the  following  tasks, 
commonly  accepted  correlations  will  be  discussed  and  comparisons  will  be  made  to 
available  data  noting  deficiencies  in  the  relations. 

When  assessing  the  various  engineering  design  codes,  the  integral  technique  ap- 
pears to  have  received  the  most  attention.  Essentially,  this  prediction  scheme  is  based 
on  the  incompressible  momentum  boundary  layer  equation  which  is  subsequently  modi- 
fied for  compressibility  by  a density  transformation  (such  as  the  Eckert  reference 
enthalpy  method).  Consider  then  the  momentum  integral  equation  for  an  irrotational 
flow,  namely 


w 


0 u 
e e 


Id,, 


f>  U 

e e 


dP  A „ 2l 
ds  L 0 e J 


(/>v) 


w 


("U), 


For  a sharp  cone  with  and  dP/ds  equal  to  zero,  one  has 
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where  the  subscript  refers  to  a zero  mass  transfer  state. 

Keeping  in  mind  that  closed-formed  solutions  are  desirable  for  expedient  para- 
metric studies,  two  concepts  are  available  for  design  purposes.  One  is  to  assume  a 
velocity  power  law  together  with  the  Blasius  form  of  skin-friction  and  a reference 
enthalpy  method  for  compressibility.  Such  a procedure  yields  the  result 


where  and  are  constants,  n and  m are  related  to  the  power  law  selected  and  the 
starred  properties  are  based  on  a reference  state  which  is  an  average  type  between 

the  wall  and  boimdary  layer  edge  conditions: 

♦ * 

p and  P = function  (h*) 
for 

h*/h  = 1/2  ( 1/2  h /h  ^ 0.  Hr  (7-  1)  M ^ 
e we  e 

The  second  approach,  which  has  been  adopted  by  many  authors,  is  somewhat  more 

sophisticated  and  attempts  to  model  more  detail  of  the  turbulent  characteristics  of  the 

boundary  layer.  If  one  consideres  either  the  von  Karman  or  Prandtl  forms  of  local 

shear  distribution  (the  Prandtl  form  will  be  adopted  here)  there  results 

- 2 2 - 
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For  a constant  pressure  boundary  layer  together  with  the  Crocco  linear  temperature- 
velocity  relation  and  the  equation  of  state,  the  density  is  related  to  the  velocity  giving 
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In  this  manner,  skin-friction  results  can  be  generated  for  the  non-blowing  condition. 

In  both  cases  presented  above,  the  design  engineer  has  resorted  to  correlations 
of  experimental  data  to  determine  the  attenuation  in  skin-friction  as  well  as  the  bound- 
ary layer  characteristic  thicknesses  ,6,  **,  0 ) resulting  from  mass  transfer.  How- 
ever, it  should  be  noted  that  the  analytical  solutions  suggest  the  form  of  the  correlations. 

In  viewing  these  attendant  difficulties  to  maintain  simplicity  in  design  codes,  the 
features  of  the  boundary  layer  characteristics  that  will  require  careful  examination 
consist  of: 

(a)  Boundary  layer  profile  shape  (velocity  power-law  exponent) 

(b)  Energy  profile  shape  (Crocco  temperature -velocity  relationship) 

(c)  Viscous  layer  thicknesses  (6 /a  , a / t,  , o/ b ) 

O 0 0 

(d)  Skin-friction  reduction  (C^C^  ) 

o 

(e)  Heat  transfer  blockage  (St/St^) 

In  particular,  the  above  will  be  examined  relative  to  improved  engineering  design 
capability  for  turbulent  boundaiy  layer  flows  that  consider  the  effects  of  angle  of  attack 
and  Reynolds  number  as  well  as  mass  transfer. 

Because  of  the  extensive  scope  of  this  study  which  examines  the  turbulent  boundary 
layer  features  listed  above,  the  body  of  this  report  is  subdivided  into  four  major  sec- 
tions. These  sections  are  listed  below  in  their  order  of  appearance. 

Section  2.  Summary  of  STREET-G  Experiment 

Section  3.  Crocco  Temperature-Velocity  Relationship 

Section  4.  Turbulent  Viscous  Layer  Thicknesses 

Section  5.  Turbulent  Heat  Blockage  and  Skin-Friction  Reduction  Due  to 
Mass  Injection 


The  STREET-G  experiment  was  summarized  because  it  represents  a unique  set 
of  data  which  heretofore  was  not  available  at  hypersonic  conditions  and,  appropriately, 
is  emphasized  in  all  phases  of  this  study. 
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SECTION  2 

SUMMARY  OF  STREET-G  EXPERIMENT 

The  experimental  investigation,  which  provides  the  main  data  base  of  the  present 
study,  represents  a two  year  effort,  the  objective  of  which  was  to  e^qperimentally 
provide  detailed  turbulent  boundary  layer  surface  and  profile  data  on  a slender  cone 
with  mass  transfer  effects.  The  study  represented  a comprehensive  investigation 
of  the  boundary-layer  mean  flow  characteristics  whose  final  product  was  a detailed 
documentation  of  the  data  including  an  assessment  of  pertinent  trends  and  observations. 
The  data  obtained  are  to  serve  two  functions  which  include:  (1)  long  term,  data  base 
provision  to  validate  sophisticated  3-D  finite  difference  or  integral  viscous  boundary 
layer  computer  codes;  and  (2)  short  term,  data  base  provision  to  validate  and  i^grade 
engineering  design  codes  which  currently  rely  upon  two-dimensional,  low  Mach 
number  data  correlations.  These  documents  are  included  as  References  1 and  2 of 
this  report. 

The  ejqjerimental  effort  concentrated  on  surface  and  boundary  layer  behavior 
relative  to  the  effects  of  mass  addition  (with  uniform  and  non-uniform  distributions), 
molecular  weight,  and  angle  of  attack  at  a fixed  Reynolds  number  condition.  The 
ii^ectant  gases  consisted  of  air,  helium,  argon,  and  Freon  with  blowing  values 
^06  = (pv)  u /(/’u)  of  .0005,  .0015  and  .0035  at  angles  of  attack  of  0®,  3®,  5®  and 

QO 

10® . The  basic  model  consisted  of  a 5®  half  angle  porous  cone  which  ^-as  60  inches 
long  and  was  constructed  from  sintered  stainless  steel  so  that  ablation  could  be  simu- 
lated by  the  injection  of  gases  through  the  porous  model  walls.  Tables  I and  n give  a 
detailed  listing  of  the  type  of  profile  data  obtained  in  the  two  year  period  while  Tables 
m and  IV  provide  a detailed  listing  of  the  surface  data. 

2.1  EXPERIMENTAL  APPARATUS 

The  test  program  was  conducted  in  the  Arnold  Engineering  Development  Center 
(AEDC),  Von  Karman  Gas  Dynamics  Facility,  Tunnel  B at  Mach  8.  This  facility  is  a 
continuous  flow,  hypersonic  wind  tunnel  with  a 50-inch  diameter  test  section.  An 
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axisymmetric  contoured  nozzle  provides  flow  at  Mach  8 v^oh  corresponds  to  free 
stream  unit  Reynolds  numbers  from  0.30  x 10®  to  3. 8 x 10®  per  foot.  The  model  used 
in  this  investigation  was  a 5 degree  half  angle  cone  which  is  nominally  60  inches  long 
(10.5  inches  base  diametex) . tt  consists  of  an  impervious  sharp  nose  (Rjij  ~ .002  is) 
which  is  9.47  inches  long,  a porous  frustum  which  is  50  inches  long,  and  a base  plate 
which  is  0.5  inches  thick  (see  Figure  1).  The  porous  section  of  the  model  is  sub> 
divided  into  fo’ic  independent  compartments,  each  with  separate  supply  tubes  which 
are  1/4  inch  0.0.  stainless  tubing.  The  forward  two  chambers  each  have  one  siq>ply 
tube;  the  two  rear  chambers  each  have  two  supply  tubes.  The  porous  frusta  of  the 
model  are  made  of  series  316  sintered  stainless  steel  (using  approximately  5 mil 
spherical  powder)  which  is  nominally  60%  dense.  The  wall  is  approximately  1/4  inch 
thick.  These  porous  elements  are  siq)ported  by  a stainless  steel  substructure  which 
contains  all  of  the  necessary  internal  plumbing  and  instrumentation  routing  ports. 

The  cone  was  instrumented  with  26  heat  transfer  gages  and  34  pressure  orifices. 
The  surface  heat  transfer  was  measured  with  a copper-constantan  Garden  type  heat 
gage.  This  gage  operates  on  the  principle  that  the  heat  flux  input  into  the  thin  sensing 
foil  is  proportional  to  the  difference  in  temperature  between  the  center  of  the  foil  and 
its  edge  (where  the  foil  edge  is  bonded  to  a copper  heat  sinl0 . Heat  transfer  measure- 
ments were  made  with  the  model  wall  at  both  ambient  temperature  and  at  an  elevated 
temperature  (the  latter  while  profile  data  were  obtained) . For  the  non-blowing  case, 
with  the  model  at  an  equilibrium  wall  temperature  (approx.  1000”  F),  the  surface  heat 
transfer  was  approximately  0.40  BTU/ft^-sec.  To  measure  this  small  value  of  heat 
transfer,  a gage  with  an  extremely  thin  foil  and  a relatively  large  diameter  was  re- 
quired. Both  of  these  characteristics  were  undesirable  for  thick  wall  porous  model 
instrumentation.  Consequently,  a gage  with  a diameter  of  1/4"  and  a constantan 
diaphragm  thickness  of  5 mils  \vas  selected.  Each  sensor  has  the  capability  of 
measuring  heat  transfer  (from  a pre-determined  calibration)  as  well  as  wall  tempera- 
ture. Surface  pressure  w'as  recorded  from  conventional  wall  static  pressure  ports. 
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Boundary  layer  profile  mean  flow  measurements  were  made  with  four  discrete 
probes;  namely»  pitot  pressure,  total  temperature,  mass  concentration  ^or  foreign 
gas  injection),  and  static  pressure  (via  a cone  probe).  Each  of  these  probes  was 
connected  to  a strut,  i^diich  in  turn,  was  connected  to  the  AEDC  probe  actuation 
mechanism.  The  probing  was  performed  in  general  with  the  pitot  and  total  tempera- 
ture probes  straddling  the  model  centerline.  The  probe  vertical  hei^t  at  y « 0 was 
ascertained  by  detecting  continuity  on  each  of  two  low  voltage  circuits  when  the  probes 
touched  the  model  surface.  Since  the  lateral  displacement  of  each  of  the  probes  (in 
the  transverse  direction)  relative  to  the  model  radius  is  small,  one  would  expect  the 
data  obtained  at  smail  values  of  the  transverse  coordinate  to  be  the  same  as  the  data 
with  each  of  the  probes  on  the  centerline  of  symmetry.  This  was  verified  ejq)erimen- 
tally  and  found  to  be  true. 


The  pitot  pressure  was  measured  with  a flattened  tip  hypodermic  tube  tip  height 
opening  approximately  0.020".  The  probe  was  designed  to  have  a small  probe  tip 
hei^t  (y)  relative  to  the  local  boundai’y  layer  thickness  (i.e.,  yA  <0.050),  to  mini- 
mize any  influence  of  the  probe  size  on  the  resulting  measurements. 

The  mass  concentration  probe  has  the  same  external  geometry  as  the  pitot  probe. 
This  probe  ;vas  connected  to  an  evacuated  gas  analyzer  system  which  contains  an 
"alphatron"  ionization  gage.  Alpha  pax’ticles  emitted  from  this  gage  ionize  the  gas 
molecules  x\’hich  are  collected  on  a cliarged  plate  and  produce  an  ionization  current. 
Data  are  obtained  by  injecting  a gas  sample  at  a prescribed  test  point  and  determining 
the  mixture  concentration.  From  pre-test  calibrations  one  can  establish,  with  this 
system,  the  mole  fraction  x \ of  the  mixture  for  each  test  point  from  the  ionization 
potential  of  the  gas  mixture. 


The  stream  total  temperature  was  measureti  with  a singly  shielded  thermocouple 
(tip  height  approximately  0.050").  This  probe  was  calibrated  in  the  free  stream  for 
several  values  of  the  Reynolds  number  and  for  a constant  total  temperature.  Additional 
calibration  points  were  obtained  from  measurements  with  the  probe  at  the  boundary 
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layer  edge  and  for  several  values  oi  the  cone  angle  of  attack.  These  data,  although 
calibrated  for  a limited  Mach  number  range,  were  assumed  to  apply  for  all  M in  the 
boundary  layer. 

One  particular  item  of  interest  in  this  overall  sfidy  u-as  to  examine  the  boundary 
layer  assumption  that  (ap/  'ly  )];j^  1 =0.  Consequently,  static  pressure  profiles  were 
measured  in  the  boundary  layer,  busmuch  as  relatively  thin  boundary  layers  would 
be  encountered  on  the  cone  (i.e.,  ^ 0.5"  for  = 0),  conventional  ogive-cylinder 
probes  would  be  unsuitable  because  cf  attendant  wall  interference  effects.  As  a 
result,  a pointed  cone  probe  was  adapted.  In  principle,  one  can  determine  the  local 
flow  properties  such  as  the  Mach  number  and  static  pressure  without  making  the 
classical  assumption  that  (i.p/dy)  across  the  boundary  layer  is  zero  (i.e.,  from  the 
measurement  of  the  pressure  on  the  surface  of  the  cone  probe  and  the  pitot  pressure 
at  that  point  and  using  conical  flo\v  theory  and  normal  sho  relations) . However,  in 
practice,  the  accurate  measurement  of  the  static  pressure  in  a thin  boundary  layer 
for  a hypersonic  flow  is  no  small  task.  Recognizing  these  difficulties,  a careful 
assessment  concerning  probe  geometry,  lateral  spacing,  and  application  relative  to 
obtaining  quantitative  data  was  made. 

The  probe  lateral  spacing  ^vas  such  that  at  a Mach  number  greater  than  approxi- 
mately 1.5,  the  cone  probe  would  not  interfere  with  the  total  temperatiure  probe  \vhich 
is  immediately  adjacent  to  it.  For  the  turbulent  boundary  layer  measurements  made 
in  Reference  1,  this  r.uitunl  interference  did  not  present  a problem.  However,  during 
the  tests  of  Reference  2,  where  tejjts  were  made  predominantly  at  a lower  Reynolds 
number  witli  measurements  made  in  the  laminar/transitional  boundary  layer  region, 
interference  effects  were  noted.  ,\s  a consequence,  the  cone  probe  was  removed  from 
the  probe  assembly  for  the  majority  oJ  the  profile  measurements. 

Readers  interested  in  more  details  of  these  probe  measurements  as  well  as  model 
and  instrumentation  details  should  consult  References  1 and  2. 


In  general,  wo  categories  of  data  were  obtained  in  this  test  program,  one  in- 
volving transient  heat  transfer  data,  and  the  second  pertaining  to  boundary  layer  pro- 
file measurements.  In  tiie  former,  the  model  was  pulsed  into  the  tunnel  flow  for  a 
period  of  approximately  five  seconds  and  heat  transfer  data  were  obtained.  The  model 
was  then  retracted  and  cooled  for  the  next  test  condition.  For  the  profile  data  entry, 
the  model  was  inserted  into  the  tunnel  flow  and  allowed  to  absorb  heat  until  an  equili- 
brium wall  temperature  was  achieved  (i.e.,  = %onv  ^ reached  thermal 

equilibrium,  mean  profile  measurements  were  obtained.  Boundary  layer  mean  flow 
profEe  data  were  obtained  for  air.  Freon- 12,  helium,  and  argon  injection  on  the  wind- 
ward and  leeward  planes  of  symmetry.  Measurements  were  made  at  several  measur- 
ing stations  at  a = 0,  ± 3“ , ± 5® , and  ± 10® . 

For  the  prescribed  free  stream  conditions  and  model  mass  injection  rates  used 
ir.  the  ejq)erimental  study,  the  pitot  pressure,  total  temperature  and  injectant  mole 
fraction  were  measured  across  the  boundary  layer  at  several  u:dal  stations.  In  addi- 
tion to  these  measurements  in  the  boundary  layer,  the  static  pressure,  heat  transfer 
and  the  temperature  at  the  model  wall  were  also  measured.  Inasmuch  as  local  flow 
properties  deduced  from  the  measured  data  are  desired,  equations  describing  param- 
eters of  a binary  gas  mixtime,  boundary  layer  characteristic  distances  and  boundary 
layer  properties  are  required.  In  particular,  the  parameters  required  are  the  mixture 
mass  fraction,  Cj,  effective  gas  constant,  R,  and  the  effective  ratio  of  specific  heats, 
7;  for  the  injectant  gases  of  helium,  Freon,  and  argon.  The  boundary  layer  param- 
eters consisted  of  the  local  Mach  number,  velocity,  static  temperature  and  density, 
viscosity  of  a binary  mixture  and  the  integral  viscous  layer  thicknesses  (a*,  o ) . 

The  equations  describing  these  parameters  are  given  in  References  1 and  2. 


11 


SECTION  3 


CllOCCO  TEMPERATURE-VELOCITY  RELATIONSHIP 


3.1  BACKGROUND 


In  recent  years  an  extensive  amount  of  work  has  been  undertaken  for  the  purpose  of 
predicting  the  behavior  of  compressible,  turbulent  boundary  layers.  This  work  has 
consisted  of  both  sophisticated  finite-difference  type  solutions  as  well  as  integral  type 
codes,  with  the  latter  demonstrating  a distinct  advantage  with  respect  to  computational 
time.  In  examining  the  basic  structure  of  the  prediction  techniques,  one  notes  a signifi- 
cant number  of  simplifications  which  are  required  for  reasonable  solutions.  Moreover, 
there  is  a very  limited  amount  of  high  speed  data  available  for  the  purpose  of 
evaluating  the  prediction  schemes. 

One  of  the  more  popular  concepts  for  predicting  compressible  turbulent  boundary 
layer  behavior  is  to  extend  the  incompressible  mixing-length  formulation  of  Prandtl  and 
von  Karman  to  compressible  flows  as  in  the  classic  woric  of  van  Driest.  Consequently, 
when  density  variations  occur  and  kinetic-heating  effects  become  appreciable,  a relation 
between  density  and  velocity  is  required.  One  such  method  of  relating  the  density  and 
velocity  is  to  assume  unit  Prandtl  and  Schmidt  numbers  throughout  the  boundary  layer. 
Often  referred  to  as  the  Reynolds  analogy,  this  assumption  'Tiplies  a relation  between 
shear  stress,  heat  transfer,  and  mass  transfer,.  The  linear  expression  relating  the 
temperature  (enthalpy)  and  velocity  fields  is  the  classic  Crocco  relation,  namely; 


H 


w 


H -h 
e w 


The  appendix  of  this  report  presents  the  pertinent  assumptions  and  a formal  derivation 
of  similitude  between  temperature  and  velocity.  Essentially,  the  equation  above  has 
been  traditionally  used  to  relate  the  density  \ ariations  in  a boundaiy  layer  to  the  velocity. 
In  so  doing,  one  must  postulate  an  equation  of  state  together  with  the  boundary  layer 
assumption  for  zero  pressure  gradient  normal  to  the  flow  direction. 
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The  development  of  the  Crocco  relationship  was  based  on  the  two-dimensional 

boundary  layer  equations  and  is  subject  to  the  constraints  of  unit  Prandtl  number,  zero 

pressure  gradient  (dP/ds),  and  constant  wall  temperature.  Hence,  the  use  of  the  Crocco 

relation  to  three-dimensional  flows  and  flows  with  mass  transfer  must  be  carefully 

examined  to  determine  its  range  of  applicability  for  engineering  prediction  techniques. 

The  Crocco  expression  will  be  evaluated  using  available  experimental  data,  emphasizing 

1 2 

in  particular  the  recently  acquired  experimental  data  of  Martellucci  and  Laganelli  ’ at 
hypersonic  conditions. 

There  have  been  several  recent  studies  concerning  experimental  verification  of  the 
Crocco  temperature-velocity  relation.  The  experiments  can  be  categorized  into 
two  types;  (1)  ’’flat  plate”  type  flows  (e.  g. , flat  plates,  sharp  cones,  and  hollow 
cylinders)  where  the  upstream  pressure  gradient  is  considered  negligible  and  (2)  nozzle 
wall  type  flows  where  an  upstream  pressure  gradient  exists.  In  the  former,  difficulty 
in  obtaining  fully  developed  turbulent  boundary-layer  flow  at  hypersonic  speeds,  and 
consequent  small  boundary  layer  thickness,  has  Impeded  experimental  studies.  On  the 
other  hand,  nozzle  walls  have  been  extensively  surveyed  since  more  accurate  measure- 
ments can  be  made  in  such  a comparatively  thicker  boundary  layer.  While  the  nozzle 
wall  features  a flow  that  simulates  expansion  over  a blunt  nose  body,  it  does  not  allow 
for  a mathematically  tractable  solution  for  flat  plate  type  flow. 

As  previously  noted,  several  experimental  investigations  have  been  performed  that 
examined  the  Crocco  velocity-temperature  relation  (References  3,  4,  and  5).  These 
works  were  critically  reviewed  by  Bushnell  et  al  who  also  performed  an  experimental  in- 
vestigation on  nozzle  walls.  The  work  of  Refex-ence  6 is  summarized  in  Figures  2 and  3. 

Here,  two  distinct  charactei'istics  are  observed  when  relating  the  data  to  the  Crocco 
variables.  In  Figure  2,  the  total  temperature-velocity  data  are  shown  for  "flat  plate” 
type  flows  (negligible  pressure  gradient)  with  significant  scatter  about  the  linear  relation. 
This  scatter  is  attributed  to  inaccurate  measurements  as  a consequence  of  small  boundary 
layer  thicknesses.  Moreover,  the  short  length  of  the  models,  compared  to  the  length  of 
wind  tunnel  nozzle  walls,  produced  transition  effects  at  the  various  measuring  stations 


which  also  contributed  to  the  data  scatter.  The  authors  note  a relatively  good  agreement 
about  the  linear,  unit  Prandtl  number  distribution.  The  data  generally  fall  in  a band  of 
.125%  of  the  linear  Crocco  distribution. 

In  Figure  3 the  Crocco  variables  are  shown  for  both  a linear  (Pr  = unity)  and  a 
quadratic  distribution  in  velocity.  Here,  the  data  appear  to  scatter  about  the  quadratic 
variation.  The  total-temperature- velocity  data  in  this  case  were  obtained  on  nozzle  wall 
type  flows  which  characteristically  display  upstream  pressure  gradient  histories.  In 
both  cases  no  apparent  effect  of  Reynolds  number  or  Mach  number  was  evident. 

Since  the  work  of  Reference  (5,  another  important  extension  of  the  Crocco  relation 
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has  I’esulted  from  the  experimental  and  analytical  investigations  of  Fiore  who  examined 
contoured  and  conical  nozzles.  Two  distinct  regions  were  observed,  neither  of  which 
followed  the  linear  Crocco  relation.  In  the  region  close  to  the  boundary  layer  edge,  the 
data  appeai'ed  to  follow  the  parabolic  velocity  distribution;  while  near  the  wall,  the 
data  tended  to  follow  a linear  distribution  with  an  off-set  slope.  The  deviation  from  the 
classic  Crocco  I'elation  was  attributed  to  a combination  upstream  history  effect  and 
longitudinal  pressure  gradient. 
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Beckworth  et  al  performed  an  experimental  investigation  on  the  wall  of  an 
axisymmetric  contoured  nozzle.  In  that  case,  it  was  noted  that  a portion  of  the  data 
tended  toward  the  cjuadratic  distribution  as  a consequence  of  the  boundary  layer  develop- 
ment from  a cold  wall  temperature  condition  at  the  throat  of  the  settling  chamber.  The 
linear  trend  of  another  portion  of  the  data  (u/u^^.  - 0.  8)  in  these  experiments  was  attributed 
to  the  thick  viscous  sublayer  that  extended  out  to  values  of  u/u^  ^ 0.  8.  It  was  also 
determined  that  the  local  density  levels  as  well  as  flow  histoiy  must  be  known  to  make 
theoretical  predictions. 
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In  a more  recent  publication,  Gates  made  a very  comprehensive  literature  review 
and  experimental  investigation  which  was  concerned  with  the  effect  of  upstream  conditions 
on  the  dov/nstream  state  of  flow.  The  experiment  employee  a half  nozzle  and  a flat  plate 
with  the  nozzle  throat  and  flat  plate  leading  edge  being  temperature  controlled.  Again  the 
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importance  of  upstream  heat  transfer,  with  and  without  longitudinal  pressure  gradient, 
was  observed.  The  energy  removed  from  the  upstream  boundary  layer  tended  to  shift 
the  downstream  enthalpy  profile  with  negligible  effect  on  the  corresponding  velocity 
p. -'file.  For  a favorable  pressure  gradient,  the  nozzle  produced  a fuller  velocity  profile 
than  the  flat  plate  case  but  relaxed  to  the  flat  plate  downstream  characteristic  conditions. 

It  was  noted  that  the  enthalpy  profile  tended  to  respond  to  local  conditions  near  the  wall 
(sublayer) , whereas  the  energy  profile  responded  to  thermal  history  effects  primarily 
in  the  turbulent  outer  region  of  the  boundary  layer.  Finally,  the  recovery  factor  for  a 
cooled  flat  plate  leading  edge  was  less  than  that  of  the  adiabatic  (nearly)  leading  edge 
case. 

Demetiiades  and  Laderman^^  have  presented  mean  profile  data  in  the  Crocco 
variables  using  the  model  and  the  same  basic  tunnel  conditions  of  References  1 and  2. 

While  the  emphasis  of  their  work  was  to  obtain  fluctuating  measurements  in  the  boundary 
layer,  mean  profile  data  were  also  obtained.  The  Crocco  energy-velocity  relation  for 
both  the  linear  and  quadratic  distributions  was  compared  with  data  including  non-blowing, 
blowing,  and  angle-of-attack  effects.  The  authors  noted  a small  effect  of  blowing  on  the 
measured  profiles.  Moreover,  at  a free-stream  unit  Reynolds  numberof  1.7X10  /ft.  the 
data  appeared  independent  of  blowing,  while  the  data  in  the  laminar  region  tended  toward 
linear  distribution  and  the  turbulent  regime  tended  toward  the  quadratic  behavior.  While 
the  authors  noted  that  the  validity  of  the  measurements  is  questionable,  there  appears 
to  be  some  irregularity  in  the  probe  data  per  se. 

As  previously  noted,  the  authors  used  the  same  model  and  tunnel  conditions  as 
were  used  in  the  experiments  of  References  1 and  2.  When  examining  the  data  of  Reference 
10,  it  is  noted  that  the  lack  of  temperature  "overshoot"  in  the  outer  regions  of  the  boundaiy 
layer  render  the  total  temperature  probe  design  questionable.  Further  comments 
concerning  the  experiments  of  Reference  10  will  be  given  in  the  text  of  this  report. 

In  a very  interesting  study,  Reda^^  conducted  an  experimental  program  that  examined 
the  effects  of  roughness  and  wavy  walls  on  turbulent  boundary  layer  skin  friction  and 
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velocity  profiles  in  a compressible,  adiabatic  flow.  Flat  and  contoured  nozzle  tunnel 
walls  were  used  in  the  experiment.  It  was  found  that  the  Crocco  relationship,  for 
adiabatic  flow,  could  be  used  to  describe  measured  temperature  profiles  for  both  smooth 
and  rough  walls. 

In  assessing  the  analytical  and  experimental  studies  of  the  Crocco  temperature- 
velocity  relation,  it  becomes  apparent  that  two  distinct  characteristics  of  the  velocity 
distribution  are  exhibited.  For  flat  plate  tj'pe  flows  (geometries  with  negligible 
longitudinal  pressure  gradients),  the  energy  distribution  tends  to  follow  a linear  velocity 
relation;  whereas,  for  tunnel  wall  (nozzle)  type  data,  the  energy  distribution  favors  a 
quadratic  type  velocity  distribution,  it  should  be  emphasized,  however,  that  the  above 
phenomena  are  not  all-encompassing.  For  example,  it  has  been  demonstrated  that  the 
upstream  thermal  history  has  some  effect  on  the  downstream  boundary  layer  characteris- 
tics. Although  the  Crocco  I’clation  does  not  appear  to  vary  with  a Reynolds  number  or 
Mach  number  variation,  the  state  of  the  boundary  lay:‘r  should  be  included  in  any 
analysis.  The  history  effect  can  be  observed  when  considering  a profile  measuring 
station  in  the  vicinity  of  the  position  of  boundary  layer  transition.  Here,  the  boundary 
layer  has  not  been  fully  developed  and  exhibits  an  upstream  history  effect  on  the  local 
properties. 

As  developed  in  the  appendix  of  this  report,  the  Crocco  linear  temperature  (enthalpy) 
velocity  r-elation  is  given  by 
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which  is  subject  to  the  constraints  of  unit  Prandtl  number,  zero  longitudinal  pressure 
gradient,  and  constant  wall  temperature.  Moreover,  the  above  relation  does  not  contain 
an  explicit  function  for  mass  transfer,  but  one  which  enters  implicitly  through  the  boundary 
conditions.  Inasmuch  as  similitude  between  the  thermal  and  velocity  fields  does 
exist  throughout  a tui'bulent  boundary  layer  (which  would  imply  that  a Reynolds  analogy 
exists)  departures  from  the  linear  lelation  are  expected.  Recognizing  that  the  data 


1« 


tended  to  deviate  from  the  linear  Crocco  relation,  Danberg**  suggested  the  following 
empirical  relation: 
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In  the  above,  for  an  adiabatic  wall  condition  0 approaches  zero  and  the  energy 
profile  is  characterized  by  a quadratic  velocity  relation  which  is  typical  of  nozzle  wall 
type  data.  However,  one  must  keep  in  mind  that  these  data  are  sensitive  to  upstream 
history  effects  (pressure  gradient  and  thermal  effects).  On  the  other  hand,  for  a unit 
Prandtl  number  (T  = T°  ),  0 approaches  unity,  and  the  energy  profile  is  characterized 
by  the  linear  distribution. 
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Kutateladze  and  Leont'ev  also  noted  that  departures  from  similarity  (which  are  a 

2 

consequence  of  the  kinetic  energy  term,  (Pr-1)  u , in  the  energy  equation  as  noted  in 

2 

the  appendix)  could  be  accounted  for  by  the  following  modification 
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The  authors  allowed  the  recovery  factor  to  be  constant,  such  that  the  parameter  ( 
was  uniquely  defined  by  Eqs.  (4)  and  (5).  Moreover,  the  parameter  5 which  is  a 
function  of  the  normal  coordinate  depends,  in  general,  on  the  longitudinal  pressure 
gradient,  the  temperature  difference  imposed,  and  the  magnitude  of  mass  transfer.  For 
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a constant  value  of  the  recovery  factor,  Eqs.  (-1)  and  (5)  can  be  expressed  in  terms  of 
the  Crocco  values,  giving 
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One  immediately  recognizes  the  value  of  the  work  of  Kutateladze  and  Leont'ev  in  the 
addition  of  the  parameter  ^ . Also,  the  above  c., nation  was  later  empirically  verified 
by  Danberg  for  the  condition  of  a unit  value  of  i,  . 


The  objective  of  the  pi-escnt  study  is  to  investigate  the  Crocco  temperature- 
velocity  relation  relative  to  hypersonic  flow  conditions.  The  examination  will  include  the 
effects  of  mass  transfer,  angle  of  attack,  and  molecular  weight  (i.  e. , with  the  injectant 
gas  differing  from  the  fx’oe  stream  gas).  Moreover,  the  study  will  emphasize  the 
STREET-G  experimental  data  together  with  data  from  the  scientific  community.  Fin- 
ally, the  suggested  temperature -velocity  relation  of  Kutateladze  and  Leont'ev  will  be 
examined  x\ith  emphasis  on  the  non-similarity  parameter  i . 


3.2  DATA  PRESENTATION 

In  order  to  assess  the  data  and  data  reduction  techniques.  Figure  1 was  constructed 
to  show  the  velocity  and  temperature  profile  data  compared  with  results  from  the 
ec|uilibrium  non-similar  boundary'  layer  program  outlined  in  Reference  13.  The  compari- 
son was  made  for  laminar  boundary  layer  conditions  with  and  without  mass  transfer. 

For  the  non-blowing  case  om-  notes  excellent  agreement  except  for  the  region  near  the 
wall  (y  < 0. 1")  whv'i'c  probe  ^suii'ace  flow  interactions  are  present.  For  the  blowing 
case  the  agreement  is  also  good  for  both  \elocity  and  temperature  profiles  except  for 
the  region  near  the  wall  (y<  0. 1")  where  once  again  inteilcixmce  effects  are  present. 

It  is  also  of  interest  to  note  that  the  measuix'd  total  temperature  profile  shows  a 5% 
overshoot  whereas  the  theoretical  value  corresponds  to  an  overshoot  of  less  than  1%. 

The  reason  for  this  diffeixmce  is  not  readily  dxvious.  It  is  suspected  that  it  is  due  to  the 
numerical  procedures  in  the  boundary  layer  program.  Nevertheless,  relatively  good 
overall  agreement  was  achieved  bi'twcen  theoiy  and  cx{X'riment. 
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When  assessing  the  available  data  in  the  literature  concerning  the  Crocco 
temperature-velocity  characteristics,  the  importance  of  upstream  history  effects  was 
noted.  In  particular,  the  relative  distance  of  the  measuring  station  to  boundary  layer 
transition  occurrence  must  be  carefully  weired.  Consequently,  Figures  5 and  6 were 
prepared  to  indicate  the  location  of  boundary  layer  transition,  with  and  without  blowing, 
for  the  zero  an^e  of  attack  conditions.  Shown  in  these  figures  are  the  axial  distributions 
of  heat  transfer  along  the  porous  cone  for  both  the  transient  and  steady  state  mode.  In 
the  latter  case,  the  model  wall  reaches  an  equilibrium  temperature  which  is  somewhat 
below  the  adiabatic  wall  value;  and  accordingly,  the  convective  heating  to  the  model 
surface  is  nominally  balanced  by  the  radiative  losses. 

There  are  two  significant  observations  that  can  be  deduced  from  Figures  5 and  6: 

first,  the  attenuation  in  heat  flux  with  blowing  and,  second,  the  effect  of  blowing  on 

transition.  In  the  former,  one  notes  a greater  heat  blockage  effect  in  the  laminar 

regime  as  compared  to  the  turbulent  case  which,  in  general,  is  compatible  with  existing 

heat  blockage  correlations.  On  the  other  hand,  the  location  of  the  transitional 

boundary  layer  regime  (defined  here  as  the  region  between  the  points  of  local  minimum 

and  maximum  heat  transfer  rates)  appears  to  be  relatively  insensitive  to  mass  transfer 

from  the  free  stream  which  is  in  direct  agreement  with  the  findings  of  Martellucci 
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et  al.  Finally,  it  should  he  noted  that  with  foreign  gas  injection,  an  increase  in  the 
heat  capacity  of  the  injectant  attenuates  the  convective  heating  level  at  the  surface. 
Moreover,  it  appears  that  the  location  of  the  transition  region  is  not  materially  affected 
by  the  frustum  injectant  gases. 

g 

Some  comparisons  of  the  transitional/turbulent  data  obtained  at  Re«  /ft  = 3.  8 x 10 

and  1.3  x 10^  for  Xoo  = 0 were  made  with  the  non-similar  boundary  layer  theory.  The 

objective  of  these  comparisons  was  to  examine  the  characteristics  of  a non-similar  theory 

with  the  validity  of  eddy  viscosity  models  for  predicting  profile  shapes.  For  convenience, 

15 

th^"*  Sraith-Cebeci  model  was  employed.  Shown  in  Figure  7 are  comparisons  of  theory 
with  velocity  and  total  temperature  profile  data  for  the  high  Reynolds  number  non-blowing 

g 

case  (i.  e. , Re«  / ft  = 3. 8 x 10  ).  The  bounds  of  the  transitional  flow  regime  as 
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deduced  from  the  surface  heat  transfer  data,  are  also  shown.  P'.xcellent  agreement 
between  theory  and  experiment  is  noted  for  this  case.  Shown  in  Figure  8 are  comparisons 
of  theory  with  the  profile  data  for  the  intermediate  Reynolds  numbers  non-blowing  case 
(Rcoo  / ft  = 1.3x10  ).  Itisinterestingtonotethepooragreementbetweentheoiyand 
experiment  for  this  case.  One  can  speculate  on  many  reasons  for  the  lack  of  agreement 
in  the  low  Rejniolds  number  condition.  Without  attempting  to  elaborate  on  these 
phenomena,  it  suffices  to  say  that  additional  work  is  required  to  better  understand  the 
applicability  of  eddy  viscosity  models  to  flows  which  are  not  "fully-developed”. 

To  illustrate  the  influence  of  mass  transfer  on  the  profile  shape,  comparisons  of 
data  for  x oo  “*  0 n ncl  A oo  — 0.  0015  (air  injection)  were  made.  Shown  in  Figure  9 are 
the  data  for  the  velocity  and  total  temperature  profiles  for  the  mass  transfer  case.  The 
non-blowing  data  are  x*epresentecl  by  the  dashed  lines  which  wei-e  faired  thru  the  data 
points.  As  was  the  case  for  the  surface  heat  ti-ansfcr,  mass  transfer  at  a constant 
blowing  rate  has  a significantly  gi'eater  influence  on  the  laminar-transitional  profile 
shape  than  for  me  turbulent  case.  Fuithermore,  it  is  of  interest  to  note  that  the  data 
obtained  at  stations  28  and  32  (i.e. , in  the  transitional  regime)  clearly  indicate  that  the 
profile  can  bo  divided  into  two  pai’t.'>.  This  observation  is  mox'e  evident  when  examining 
the  temperature  profile  data  of  Figure  U.  it  appears  that  the  laminar  boundary  layer  has 
been  "lifted"  from  the  surface  the  injei-ted  gases.  As  a result  one  obtains  a region 
near  the  wall  which  has  a quasi  constant  total  enth.ilpy  near  the  injectant  value,  while 
the  outer  part  of  the  profile  resembles  a laminar-like  profile. 

One  method  of  establishing  whethei  the  data  obtaiiu  d exliibit  a fully  developed 
characteristic  is  through  examination  of  the  \elocity  profile  exponent,  n.  The  velocity 
power  law  exponent,  n,  has  been  plotted  :.s  a function  of  Re^  for  the  zero  blowing,  zero 
angle  of  attack  case  and  is  shown  in  1- igare  it).  Fach  of  the  data  points  is  identified  as 
to  the  type  of  flow  experienced,  ue. , laminar,  transitional,  or  tuiioulent.  When  the 
data  of  Figure  10  are  compared  with  the  data  band  of  References  K;  and  17  it  is  seen  that 
the  turbulent  flow  points  (10^  ■ Re„  ■ 10*)  all  lie  in  the  "overshoot"  region  where  the 
turbulent  boundary  layer  has  not  yet  come  to  full  ecjuilibhum  (i.e. , fully  developed). 
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This  is  a rather  surprising  result  when  considering  the  model  length  of  60  inches  and 

g 

maximum  facility  Reynolds  number  of  3. 8 x 10  /ft  which  yields  over  3 feet  of 
"turbulent"  flow  (i.  e. , the  end  of  transition  as  depicted  from  surface  heat  transfer 
data  is  at  23  inches).  This  observation  raises  some  speculation  concerning  published 
"turbulent"  boundary  layer  data.  This  leads  one  to  believe  that  the  only  available 
hypersonic  "fully  developed"  turbulent  boundaiy  layer  data  may  be  those  works  obtained 
on  nozzle  walls. 

Another  item  discussed  in  the  assessment  of  the  literature  review  is  boundary  layer 
thickness.  This  becomes  of  concern  when  the  probe  size  relative  to  the  boundary  layer 
height  and  the  potential  interference  effects  of  the  probe  and  wall  may  produce  question- 
able data.  Shown  in  Figure  11  are  the  measured  hydrodynamic  thicknesses  for  several 
blowing  values  at  both  the  free  stream  unit  Reynolds  numbers  tested  in  the  program.  For 
convenience,  the  transition  zones  have  also  been  shown.  One  notes  the  significant 
fliickening  of  the  boundary  layer  with  mass  transfer  which  allows  an  extended  profile 
survey.  However,  the  total  pressure  and  total  temperature  probes  used  in  the  experi- 
ment (approximately  0.050  inches  in  diameter)  are  subject  to  erratic  readings  below 
station  35  (non-injectant  cases)  for  the  high  Reynolds  number  case  and  below  station  45 
for  the  low  Reynolds  number  conditions.  It  is  also  interesting  to  observe  the  change 
in  slopes  between  the  two  Reynolds  number  cases.  For  the  highest  condition,  the  cun’es 
are  convex  and  grow  according  to  the  4/5  power  of  distance  as  expected.  However,  the 
low  Reynolds  number  case  exhibits  curves  that  are  concave.  This  is  believed  to  be  a 
consequence  of  the  effects  of  transition  occurring  in  the  bouno,  ry  layer  over  a larger 
section  of  the  model. 

Figure  12  exhibits  the  surface  pressure  variation  as  functions  of  axial  distance.  The 
data  represent  both  uniform  and  non-uniform  blowing  conditions.  Concerning  the  latter 
condition,  one  obvious  region  of  discontinuity  occurs  at  the  juncture  of  the  impermeable 
nose  with  the  first  porous  chamber.  In  the  tests  of  Reference  2,  combinations  of  injec- 
tion through  the  individually  controlled  chambers  were  devised  to  investigate  a differential 
mass  transfer  process.  In  particular,  chambers  1 and  2 were  set  at  one  blowing 
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distribution,  while  chambers  3 and  4 were  controlled  at  twice  the  value  and  no  injectant, 

respectively.  Clearly  a history  effect  is  expected,  and  the  relaxation  phenomenon  was 

evident.  A detailed  discussion  concerning  the  non-uniform  injection  test  was  given  by 

18 

Laganelli  and  Martellucci  and  will  be  discussed  further  herein.  Concerning  Figure  12 
two  distinctive  characteristics  are  interesting  to  note  here.  First  the  step-down  in 
blowing  case  (15/15/0/0)  showed  no  effect  on  pressure  compared  to  the  step-up  case 
(15/15/30/30).  Secondly,  an  incr-ease  in  pressure  is  apparent  at  the  juncture  of  the 
impermeable  tip  and  chamber  1 as  a result  of  blowing  (s  a:  10  inches).  The  same 
effect  is  not  observed  at  the  juncture  of  discontinuity  for  the  step-up  injection  case,  but 
occurs  in  a distributed  fashion  downstream.  This  can  be  attributed  to  the  thicker 
boundaiy  layer  which  is  more  evident  at  the  chamber  2/3  juncture  than  at  the  nose 
section. 

3.3  CRQCCO  TEMPERATURE-VELOCITY  DATA  RESULTS 

The  Crocco  tempo ratuiv- velocity  data  will  be  presented  in  this  section  under  five 
conditions  of  flow.  These  include  uniform  injection  at  zero,  leeward,  and  windward 
angles  of  attack;  non-uniform  injection;  and  molecular  weight  behavior.  Interpretation 
of  these  data  will  be  provided  in  the  penultimate  section  of  this  I'epoii:.  Finally,  it 
should  be  noted  that  no  attempt  has  been  made  to  eliminate  data  points  that  are  obviously 
in  error  as  a consequence  of  probe  interference. 

3.3.1  Uniform  Injection;  a = 0° 

Figui’cs  13  and  14  represent  the  Crocco  variables  with  and  without  mass  transfer 
effects  at  the  high  Reynolds  number  conditions  of  these  tests  (Re^  /ft  = 3.  7 x 10  ). 

In  Figure  13,  several  profile  stations  are  represented  for  zero  mass  transfer  conditions. 
Here,  one  notes  that  the  data  tend  to  significantly  deviate  about  the  linear  Crocco 
relation.  For  this  condition  (i.e. , Re«,/ft  = 3.  7 x 10  ) boundary  layer  transition  onset 
and  end  are  approximately  located  between  stations  13  and  22.  Hence,  the  relative  size 
of  the  boundary  layer  thickness  together  with  transition  I'ender  the  results  questionable 
in  this  region.  On  the  other  hand,  the  data  at  stations  39  and  52,  although  not  fully 
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developed,  indicate  a quadratic  effect  prior  to  the  energy  "overshoot"  region.  This  latter 
condition  is  quite  typical  of  supersonicAypersonic  flow  behavior. 

A significant  increase  in  boundary  layer  thickness  was  noted  with  increasing  mass 
transfer  which  allows  for  a more  effective  boundary  layer  survey.  Figure  14  was 
constructed  to  observe  the  effect  of  mass  transfer  at  a fixed  station.  Here,  the  data 
tend  to  follow  a modified  quadratic  relation  with  temperature  overshoot  at  approximately 
90%  of  u/u^.  It  should  be  noted  that  for  this  figure  only,  data  points  have  been  eliminated 
that  reflect  possible  probe  interference  effects. 

Figures  15  and  16  express  the  enthalpy  (temperature)  - velocity  characteristics 

0 

for  a free  stream  Reynolds  number  of  1. 3 x 10  per  foot.  For  this  condition,  the 

transition  zone  (based  on  healing)  was  in  the  region  of  stations  20  to  30  ( Xoo  =0)  and 

stations  20  to  34  ( x «>  f 0).  In  Figure  15  one  notes  that  the  data  tend  toward  the  classic 

Crocco  linear  distribution  and  are  further  characterized  by  the  temperature  overshoot 

at  the  outer  region  of  the  boundary  layer.  Moreover,  there  appears  to  be  less  scatter  in 

6 

the  data  than  for  the  3. 7 x 10  Reynolds  number  tests  (Fig.  13)  even  though  the  corres- 
ponding boundary  layer  thicknesses  are  smaller.  In  Figure  16,  the  data  tend  to  exhibit  a 
peculiar  behavior.  In  particular,  stations  28  and  32,  which  are  the  transition  zone, 
deviate  significantly  close  to  the  wall  and  then  tend  toward  the  quadratic  distribution.  On 
the  other  hand,  the  data  in  the  downstream  region  of  the  transition  zone  appear  to  follow 
some  type  of  modified  quadratic  distribution  as  was  evident  in  the  higher  Reynolds  number 
cases.  It  is  interesting  to  note  that  the  data  in  the  transition  zone  for  the  non-injectant 
case  (Figure  15)  did  not  exhibit  the  odd  behavior  of  the  mass  transfer/transition 
coupled  case. 

Figure  17  represents  the  enthalpy-velocity  characteristics  for  a free  stream 

0 

Reynolds  number  of  0,4  x 10^  per  foot.  This  condition  allows  for  an  examination  of  data 
in  a laminar  boundary  layer  state.  0\\’ing  to  probe  Interference  effects  close  to  the  wall, 
the  data  follow  the  linear  distribution  with  the  usual  temperature  overshoot.  The  same 
pattern  is  exhibited  at  the  low  Reynolds  number  condition  as  was  experienced  at  the  higher 
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Reynolds  number  values;  namely,  mass  transfer  tends  to  alter  the  data  from  the  linear 
toward  the  quadratic  distribution. 

Figures  18  and  19  were  constnacted  to  examine  a Reynolds  number  effect  with  and 
W'it’  out  mass  transfer.  In  Figure  is,  the  data  are  shown  for  the  region  around  station  i 

j 

39.  Here,  no  discernible  effect  of  Reynolds  number  is  apparent.  If  one  is  allowed  some  ; 

license  to  speculate,  it  is  possible  that  the  laminar  flow  results  (low  Reynolds  number)  ' 

are  showing  a distinct  variation  from  the  unit  Prandtl  number  result  when  compared  with  1 

the  turbulent  and  semi-tui-bulent  data.  (The  latter,  which  is  characterized  by  the  ’■ 

*■ 

mid-range  Reynolds  number,  is  subject  to  upstream  transition  effects. ) P'or  example, 

1 

the  laminar  recove ly  factor  is  less  than  the  turbulent  value;  hence,  the  data  should 
reflect  a larger  variation  foi  a non-unit  Prandtl  number  condition.  In  Figure  19, 

the  effect  of  Reynolds  number  with  mass  transfer  is  investigated  in  the  vicinity  of  | 

station  30.  While  the  data  tend  to  deviate  fiom  the  linear  distribution,  there  does  not 
appear  to  be  any  dependence  on  Reynolds  numoer. 

i 

r 

3.3.2  Uniform  Injection  - Windward  Angle  of  .-Utack  = 

Figures  20  and  21  present  typical  v.’roeeo  <'nthalpy-vrlocity  data  for  windward  angles  ] 

J 

of  attack  at  the  high  free  stream  ReynoUls  number  condition.  An  inspection  of  these 

data  indicates  that  the  linear  relation  tends  to  fit  the  results  with  and  without  mass  i 

transfer  effects.  Th<>  data  also  displaj  the  energy  o\ershoot  characteristics;  however, 
the  level  of  overshoot  is  not  as  pi-onouiu-ed  as  for  tin  zero  angle  of  attack  case.  Moreover, 

5 there  doe-s  not  appear  to  In*  any  di'.tiiK  tion  in  the  tiata  trend  from  the  linear  distribution 

with  increasing  angle  of  attack.  Another  interesting  asp  -ct  of  these  data  is  that  no 
significant  variation  between  the  reported  stations  is  ap)parent.  This  phenomenon, 
which  was  not  characteristic  for  the  zero  angle  of  attack  data,  occurred  with  and  without 
mass  transfer  (see,  for  example.  Figure  21). 

3.3.3  Unifonn  Injection  - Leeward  Angle  of  Attack 

Typical  Icewaitl  angle  of  attack  data  ar<*  presented  in  Figui-es  22  and  23  with  and 
without  mass  transfer  effects.  A Mgnificant  departuri"  from  the  windwaixl  or  zero  angle 


of  attack  results  is  observed.  In  the  present  case,  the  data  tend  to  follow  the  quadratic 
behavior  as  well  as  exhibit  variations  between  the  measuring  stations.  Moreover,  the 
data  (Station  39)  also  show  a mass  transfer  effect,  whereby  increasing  the  mass  transfer 
parameter  \oo  shows  departures  from  the  quadratic  relation  ( x <»  = 0)  to  some 
modified  form  of  the  quadratic  relation.  Once  again,  the  level  of  energy  overshoot  is 
not  as  pronounced  as  in  the  zero  angle  of  attack  eases.  The  variation  in  data  between 
stations  39  and  52  with  (Figure  22)  and  without  (Figure  23)  blowing  is  not  known. 

3.3.4  Non-Uniform  Injection:  a =0° 

One  of  the  more  interesting  tests  performed  during  the  experimental  program  of 
References  (1)  and  (2)  consisted  of  non-uniform  mass  transfer  along  the  sharp  porous 
cone.  Since  each  of  the  four  compartments  of  the  cone  was  individually  controlled  for 
blowing  characteristics,  various  combinations  of  injection  could  be  examined.  In  this 
manner,  one  can  investigate  upstream  history  effects  primarily  due  to  mass  transfer. 
Figures  24  and  25  relate  the  Crocco  variables  to  non-uniform  injection  data.  Surface 
and  profile  data  concerning  the  same,  as  well  as  several  other  combinations  of  injection 
can  be  found  in  References  (1),  (2),  and  (18).  Figure  24  shows  the  injection  combination 
for  Xoo  =0.0015  in  chambers  1 and  2 and  Ao»  = 0 in  chambers  3 and  4.  Also  shown 
for  comparison  is  the  data  band  for  the  uniform  injection  case  of  A «=  =0.  0015. 

The  data  for  the  non-uniform  case  tend  to  follow  the  quadratic  velocity  relation  and 
proceed  through  the  energy  overshoot  region  in  the  usual  fashion.  It  is  interesting  to 
note  that  the  data  deviate  slightly  from  the  uniform  injection  case  and,  more  importantly, 
show  significant  departures  fi’om  the  zero  injection  case  (Figure  13).  The  latter  condi- 
tion indicates  that  the  boundary  layer,  relative  to  the  Crocco  variables,  does  not  tend  to 
recover  to  the  non-injection  distribution.  An  examination  of  Figure  12  shows  that  the 
pressure  distribution  recovers  to  the  zero  injection  value  thereby  indicating  some  other 
mechanism  for  the  departures.  It  should  also  be  noted  that  both  uniform  and  non-uniform 
cases  are  subject  to  turbulent  boundary  layer  conditions. 


In  Figure  25,  a non-uniform,  step-up  injection  case  is  pi*esented.  As  in  the  previous 
figure,  a data  band  for  the  uniform  injection  distribution  is  shown  for  comparison.  Here, 
one  notes  still  a different  pattern;  namely,  a departure  from  the  uniform  case  toward  the 
linear  Crocco  distribution.  These  data  are  not  significantly  different  from  the  uniform 
injection  case  (Figure  14)  for  the  stepped-up  injection  rate.  Furthermore,  Figure  12 
shows  an  increase  in  pi’essure  as  a consequence  of  the  step  in  injection  and  recovers  to 
the  uniform  lower  distribution  case.  Hence,  an  apparent  upstream  pressure  gradient 
effect  tended  to  drive  the  data  toward  the  linear  distribution  (the  opposite  of  what  may 
be  expected)  which  indicates  the  blowing  level  is  potentially  in  control  of  the  profile 
characteristics. 

3.3.  3 Molecular  Weight  Effect 

Figures  2(5  and  27  present  the  Croeeo  enthalpy- velocity  relation  for  the  situation 
where  the  injectant  gas  differs  from  the  boundary  layer  fluid.  In  particular,  injectants 
of  helium,  Freon-12,  and  argon  wei\-  examined,  'j'he  enthalpy  was  determined  from  a 
value  of  specific  heat  of  the  injectant-bouiKlary  layer  gas  mixture,  which  in  turn,  was 
calculated  from  mole  fraction  measurements.  The  results  are  difficult  to  interpret, 
but  will  be  presented  with  further  comments  in  the  next  section.  In  Figure  2ti,  a signifi- 
cant depaiture  from  the  Um  ar  relation  is  obseixed  when  the  injectant  gas  differs  from 
the  main  stream,  'I’here  is  no  trend  m the  ilepaiture  from  the  air/air  result  with  either 
molecular  weight  or  specific  heat.  Wink  the  argon  data  are  close  to  the  air  results 
(their  being  similar)  a substantial  deviation  is  observeii  in  the  helium  data, 

Freon-12  data  are  compai’cd  to  the  air  results  in  l igiire  27  for  stations  39  and  52. 
While  no  significant  departures  ai  e recognised  lx  tween  stations,  the  Freon-12  data  appear 
to  follow  the  lineal  relation  more  closely  than  the  air.  Moreover,  the  enei’gy  overshoot 
region  appears  to  be  less  magnifietl  with  the  non-air  injectant  gases. 

3.4  DISCUSSION  OF  RESUI.TS 


As  noted  in  tJie  background  sceti  i;.  o|  this  report,  sliarp  cone  t\  pe  data  tend  to 
follow  the  linear  Crocco  relation.  'I'his  i.,  a conse(|Uence  of  a negligible  longitudinal 


pressure  gradient  and  a Prandtl  number  (for  an  air  free  stream)  close  to  unity.  In 

general,  this  behavior  prevailed  in  the  tests  of  References  1 and  2 owing  to  the  devi 

in  data  as  a result  of  the  small  boundary  layer  thicknesses  as  well  as  attendant 
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upstream  transition  effects.  Moreover,  Danberg  and  Kutateladze  and  Leont'ev 
recognized  that  departures  from  the  linear  relation  could  be  attributed  to  the  kineti 
energy  of  the  stream  and  suggested  modifications  to  the  classic  linear  relation,  T1 
modified  Crocco  relation,  which  is  discussed  in  the  appendix,  is  given  by 

' e ' 

where  fl  = {h  - h )/(H  - h ).  Figure  28  shows  the  parameter  fl  as  a functio 
aw  w e w 

the  dimensionless  injection  term  X«  for  zero  angle  of  attack  conditions.  One  wil 
that  stations  39  and  52  represent  the  only  stations  for  which  data  were  reported  at 
three  Reynolds  numbers.  The  curve  was  not  drawn  to  express  a functional  form  ( 
with  Xoo  , but  to  indicate  the  range  of  kinetic  energy  parameter  for  these  tests, 
notes  that  blowing  has  an  effect  on  the  ^ parameter  which  is  expected  inasmuch 
wall  enthalpy  is  reduced  with  mass  transfer.  Essentially,  varied  between  . f 
. 75  for  the  major  portion  of  the  tests. 

Shown  in  Figure  29  are  t vo  groups  typical  of  the  data  presented  which  are  cr 
with  curves  calculated  using  E<),  (G)  for  {;  equal  to  unity  and  the  indicated  vali 
P.  Here,  one  notes  that  the  modified  quadratic  relation  is  in  concert  with  the  dat 
but  appears  to  overpredict  the  results.  Moreover,  the  modified  quadratic  relat 
with  5 equal  to  unity  does  not  compare  with  the  data  trend  in  the  energy  overs 
region.  Consequently,  Figure  30  was  constructed  using  Eqs,  (4)  and  (5)  with  tl 
for  several  cases.  Here,  the  parameter  ^ is  shown  as  a function  of  dimensic 
distance  through  the  boundary  layer.  Two  distinct  regions  are  evident:  first,  t 
region  dost  to  the  wall  where  departures  from  unity  ex<st  as  well  as  reversal  i 
presence  of  blowing;  and  secondly,  the  energy  overshoot  region  (y/6  > . 55  - . 
While  several  of  the  data  points  in  the  wall  region  are  subject  to  error,  the  dep 
from  unity  as  well  as  the  trend  with  injection  are  considered  accurate. 
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As  noted  above,  the  similarity  parameter  shows  a decrease  with  blowing  in  the 
region  close  to  the  wall.  Other  interesting  effects  on  the  similarity  parameter  ^ include 
molecular  weight  (Figure  30B),  at  fixed  values  of  X«  and  Re  , where  no  apparent 

oc 

trend  is  obvious.  Moreover,  for  X«>  eoual  to  zero,  there  does  not  appear  t ) be  any  angle 
of  attack  effect  (Figure  30C);  however,  with  blowing  the  leeside  data  indicate  a more 
significant  variation  than  the  windward  data.  Finally,  a Reynolds  number  effect  (Figure 
30D)  shows  a significant  departure  between  the  laminar  boundary  layer  and  the  combined 
transitional/turbulent  boundary  layer  results. 

With  the  exception  of  the  region  close  to  the  wall,  the  parameter  i.  appears  to 

approach  a value  equal  to  the  recipi’ocal  of  the  recovery  factor  (no  mass  transfer). 

On  the  other  hand,  the  I’ecovery  factor  is  expected  to  exhibit  more  of  an  effect  toward  the 

wall,  with  blowing,  than  at  the  edge  of  the  boundary  layer.  Hence,  for  mass  transfer  the 

pai'ameter  ^ was  assumed  to  take  the  form  f,  2; ~ /•iV  These  results  are  shown 
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in  Figure  31  whore  /3  was  assumed  to  be  1/  2,  which  is  typical  for  zero  injection  data. 
Hero,  one  notes  that  the  trend  correctly  pi'e<liets  what  is  observed  in  the  data;  that  is,  the 
modified  quadratic  relation  tends  towani  the  linear  distribution  in  tlie  presence  of  mass 
transfer.  This  trend  is  by  no  means  a necessary  consequence  of  mass  transfer;  but 
rather,  a result  of  the  temperature  potential  between  the  adiabatic  and  wall  values. 

Finally,  the  inclusion  of  the  parameter  correctly  predicts  the  energy  overshoot  behavior 
which  is  characteristic  of  high  speed  flow. 

In  the  preceding  discussion,  it  was  noted  that  the  linear  Crocco  relation  could  be 
expressed  as 


where  1 - ‘ ' 2.  Moreover,  it  was  noted  in  Reference  19  that  the  use  of  the  exponent, 

f , is  suitable  for  mass  transfi'r  conditions  as  well  as  non-injection  conditions.  A re- 
examination of  the  Crocco  relation  has  shovn  that  the  assumed  exponent  dependence  on 


mass  transfer  (which  is  actually  a boundary  condition)  was  erroneous.  This  can  be  seen 

by  an  inspection  of  Eq.  (6).  Here,  for  the  adiabatic  condition  (h  -*  h ),  fl  -♦  0 
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and  results  in  the  quadratic  case.  On  the  other  hand,  for  h _ H , fl  ui.itv, 
the  linear  case.  Inasmuch  as  a mass  transfer  process  tends  to  cool  the  wall  (hence, 
drive  toward  an  isothermal  condition,  0 1),  the  corresponding  effect  is  to 

drive  the  enthalpy  potential  toward  the  classic  linear  distribution.  The  attempt  in  this 
study  was  to  model  the  effect  of  mass  transfer  in  the  modified  Crocco  relation  throu^ 
the  non-similarity  parameter  { . 

3.5  CONCLUSIONS 

The  classic  Crocco  temperature-velocity  relation  has  been  evaluated  experimentally 
at  hypersonic  flow  conditions  on  a sharp,  porous  cone.  The  evaluation  was  made  with  the 
effects  of  mass  transfer,  angle-of-attack,  and  molecular  weight  variations  of  the 
injectant  gas  over  a range  of  free  stream  Reynolds  numbers  of  0.4  to  3.  7 x 10  /ft  and 
fixed  free  stream  Mach  number  of  approximately  8. 0.  Departures  from  the  linear  Crocco 
relation  (which  is  constrained  to  unit  Prandtl  number,  zero  longitudinal  pressure  gradient, 
and  constant  wall  temperature)  were  obsen^ed  that  tended  to  follow  a modified  quadratic 
velocity  behavior.  The  data  were  compared  to  the  suggested  modified  quadratic  relation 
of  Kutateladze  and  Leont'ev  and  appeared  to  agree  favorably,  including,  as  well,  the 
enei^y  overshoot  region. 

Several  of  the  pertinent  obseiwations  made  in  this  study  are  as  follows: 

1.  In  general,  the  data  tended  to  follow  the  linear  Crocco  relation  with  a temperature 
overshoot  resulting  from  the  kinetic  energy  of  the  flow.  This  result  is  typical 

of  flat  plate  type  flows  (sharp  cone)  which  have  previously  been  investigated. 

2.  As  a consequence  of  the  small  boundaiy  layer  thicknesses,  the  probe  data  in 
the  region  of  y <.  0. 1"  are  considered  questionable. 
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3.  A comparison  of  the  pi-esent  test  results  with  other  data  ’ ‘ substantiates  that 
the  Mach  number  does  not  appear  to  have  any  effect  on  the  temperature- velocity 
relations. 
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4.  For  these  cases  tested,  as  well  as  comparison  to  other  reported  data®’  ^ 

Rey  nolds  number  had  no  apparent  effect  on  the  Croceo  relation. 

5.  Mass  transfer  effects  (air-injection)  tended  to  depart  from  the  linear  relation 

toward  a modified  quadratic  velocity  relation.  It  is  believed  that  such  departures 

are  a consequence  of  the  increased  temperatuix)  potential  (i.e. , T - T ) as 
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well  as  changes  in  the  velocity  profile.  Essentially,  increased  blowing  lowers 
the  wall  temperature  value,  with  ii  increasing  as  well,  and  tends  toward  the 
quadratic  Crocco  relation. 

6.  As  noted  by  several  experimenters,  upstream  histoiy  effects  should  be  considered 
when  evaluating  the  Crocco  i*elation.  In  flat  plate  type  flows  the  relative 
distance  of  the  measuring  station  to  the  boundary  layer  transition  region  must 

be  considered.  In  the  pi’osent  tests,  departures  from  the  reporting  stations 
were  noted  which  can  be  attributed  to  flow  development.  The  reporting  stations 
in  the  transition  region  (Re»  / ft  = 1.3  X 10®)  per  se  indicated  interesting 
results.  Without  mass  transfer,  the  data  tended  toward  the  linear  relationship, 
while  with  blowing  effects  u significant  departure  was  observed  that  eventually 
became  quadratic. 

7.  Laderman^^^  has  reported  that  the  laminar  boundary  layer  flow  tended  to  agree 

with  the  linear  relation,  while  the  turbulent  flow  state  followed  the  quadratic 
relation.  The  data  of  the  present  study  are  too  inconclusive  to  separate  the 
state  of  flow  relative  to  the  Ci'occo  relation.  At  the  high  Reynolds  number 
conditions,  the  reporting  stations  (for  a = X fic  “ 0)  tended  toward  the 
linear  relation  except  for  station  52  (end  of  the  model)  which  displayed 
t'uadratic  characteristics.  Un  the  other  hand,  the  data  for  the  mid-range 
Reynolds  number,  for  « = x ..  =.  o,  appeared  to  agree  quite  well  with  the 

linear  relationship  which  also  characterized  the  low  Reynolds  number  tests. 

t).  Angle  of  attack  effects  on  the  Crocco  relation  indicated  a specific  pattern; 
namely,  windwai’i'.  orientation  tended  towaixl  the  linear  relation  while  the  lee- 
ward side  was  characterized  by  the  quadratic  behavior.  It  should  bo  noted 
that  the  aliove  is  completely  opposed  to  the  obseiwations  of  Reference  10. 

Keeping  in  mind  tliat  tlic  same  model  and  facility  were  used  at  slightly 
different  Reynolds  numlier,  the  difference  in  the  data  has  not  been  clarified. 
Windward  and  leeward  noincncljture  was  considei’cd  identical  based  on  the 
relative  sizes  of  tlu-  boundary  layer  thicknesses  for  the  tests  of  Reference 
10  and  tlie  present  experiments. 


' 9.  For  non-uniform  injection,  both  a step-up  and  step-down  in  blow'ing  followed  a 

' modifieil  quadratic  distriliution.  In  thi>  latter,  the  data  were  closer  to  the 

^ quadratic  relation  than  tlie  step-up  case  which  exhibited  more  nearly  linear 

I characteristics. 
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10.  No  specific  conclusions  could  be  made  for  the  non-air  injection  tests.  The 
heavy  molecular  weight  gases,  argon  and  Freon-12,  followed  the  linear 
Crocco  relation,  while  helium  injection  showed  a significant  departure 

in  the  opposite  direction  of  the  quadratic  relation.  Further  woric  is  required 
before  non-air  injection  cases  can  be  used  readily  in  Crocco  coordinates. 

11.  It  appears  that  the  modified  quadratic  Crocco  relation  suggested  by  Kutateladze 
and  Leont'ev  can  be  used  for  engineering  prediction  techniques.  In  this 
situation,  the  non-similarity  parameter  ? can  be  replaced  by  the  reciprocal 
of  the  recovery  factor;  i.  e.  i = Pr~^/3  is  recommended  for  use  in  Eqs. 

(6)  or  (7). 
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SECTION  4 

TURBULENT  VISCOUS  LAYER  THICKNESSES 
4.1  BACKGROUND 

There  are  many  problems  associated  with  hypersonic  flight  for  which  the  predic- 
tion of  characteristic  thicknesses  of  the  boundary  layer  would  be  desirable.  These 
inclULe  entropy  gradient  effects  on  convective  heating  rates,  signal  attenuation  in 
plasma  and  wake  studies,  boundarj'  layer  separation  (surface  discontinuities  and  con- 
trol surfaces),  rocket  engine  nozzle  flow,  and  acoustic  phenomena.  Information  con- 
cerning the  growth  of  the  viscous  layer  thicknesses  further  enhances  the  design  engi- 
neers capabilities  to  examine  three-dimensioruil  flows  for  complicated  shapes.  While 
our  understanding  and  analysis  of  three-dimensional  flows  has  steadily  advanced  in 
recent  years,  the  difficulties  involved  in  the  theor)'  of  turbulent  shear  flow  have  resul- 
ted in  inadequate  engineering  prediction  capability.  Moreover,  integral  techniques, 
which  generally  employ  streamline  coordinates  with  small  cross  flow  and  local  simi- 
larity. have  not  been  successfully  developed  for  turbulent  boundary  layer  flows. 

As  a consequence  of  those  attendant  difficulties,  tiie  design  engineer  has  resorted 
to  semi-empirical  correlations  based  on  available  experimental  data.  However,  these 
attempts  have  been  impaired  because  of  the  lack  of  detailed  data  as  well  as  an  incom- 
plete experimental  definition  of  the  mean  profile  ciuracteristics  across  the  boundary 
layers  of  hypersonic  flows.  Several  engineering  analyses  have  been  developed  within 
the  scientific  community  which  are  concerned  w itli  computing  the  vehicle  loads  and 
drag  (which  implies  an  accurate  definition  of  the  local  properties).  However,  an  exam- 
ination of  these  techniques  is  necessary  inasmuch  as  recently  acquired  data  has  ind'  ■ 
cated  thiit  modifications  are  recjuired  in  tlie  concepts.  In  |)articular,  the  effects 
mass  transfer  (including  non-air  injectants)  and  angle  of  attack  on  the  viscous  layer 
thicknesses  should  be  carefully  assessed. 


Inherent  with  the  above  phenomena  is  the  requirement  of  the  state  of  the  boundary 
layer  itself.  Here,  fully-developed  turbulent  flow  regions  should  be  considered  sep- 
arately from  transitional  flow  regimes  inasmuch  as  the  semi-empirical  correlations 

are  based  on  fully-developed  turbulent  flow.  Turbulent  boundary  layer  velocity  pro- 
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files  are  customarily  characterized  by  a power  law  relation  u/u^  = (y/s  ) where  the 
exponent,  n,  assumes  the  value  of  7 for  a fully  developed  turbulent  boundary  layer 
flow.  However,  an  examination  of  the  overshoot  region,  defined  as  the  end  of  transi- 
tion (where  the  point  of  maximum  heat  flux  occurs),  indicates  that  the  parameter  n is 
augmented  to  values  as  high  as  16.  Hence,  the  impact  on  the  viscous  layer  thicknesses 
resulting  from  boundary  layer  development  will  also  require  cax'eful  assessment. 

Inasmuch  as  the  emphasis  of  this  study  is  to  verify/modify  existing  semi-empirical 
techniques  based  on  experimental  data,  one  is  restricted  to  data  availability.  Recog- 
nizing this  difficulty,  Martellucci  and  Laganelli^’  ^ had  undertaken  an  experimental 
study  to  provide  detailed  turbulent  boundary  layer  profile  data  for  an  axisymmetric  con- 
figuration at  both  zero  and  non-zero  angle  of  attack.  The  model  used  for  this  investiga- 
tion was  a 5-degree  half  angle  porous  cone  which  was  60  inches  long.  Gases  consisting 
of  helium,  argon,  air,  and  Freon-12  were  used  to  simulate  ablative  conditions  at  sev- 
eral Reynolds  number  conditions.  These  experiments  were  intended  to  augment  the 
earlier  work  of  Martellucci  and  Rie^^’  which  found  a significant  departure  from  ex- 
isting industrial  code  prediction  capability  to  their  experimental  data. 

An  example  of  this  disparity  is  depicted  in  Figures  32  through  34.  It  should  be 
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noted  that  the  data  of  Reference  20  (open  symbols)  and  the  flat  plate  data  of  Danberg' 
(filled  symbols)  are  shown  together  with  zero  angle  of  attack  industrial  prediction 
techniques.  It  is  evident  that  each  of  the  methods  overpredicts  the  effect  of  mass 
transfer.  Moreover,  contrary  to  the  assumptions  made  in  these  prediction  techniques, 
the  form  factor  h /b  does  vary  with  mass  transfer  rate.  Although  several  early  the- 
ories  for  turbulent  boundary  layer  with  mass  addition  were  available,  ’ ’ expei’i- 

mental  skin  friction  data  did  not  agp'ee  with  theories.  ^ Of  paramount  interest, 
however,  is  the  fact  that  the  industrial  codes  relied  heavily  on  these  early  works  and 
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often  did  not  change  as  new  information  became  available.  Consequently,  the  early 
theories  are  not  necessarily  considered  to  provide  adequate  descriptions  for  turbulent 
boundary  layer  behavior. 

The  prediction  of  the  turbulent  viscous  layer  thicknesses  at  angle  of  attack  relies 
on  empirical  curve  fits  of  limited  data  (in  the  absence  of  mass  transfer).  More- 

over, the  additional  complication  of  mass  transfer  coupled  with  angle  of  attack  on  pre- 
diction capability  is  virtually  non-existent. 

As  previously  noted,  the  state  of  the  boundary  layer  relative  to  its  fully-developed 
characteristics  should  be  known  prior  to  assessing  viscous  layer  thicknesses.  A sur- 
vey^® of  velocity  profile  data  from  turbulent  boundaiy  layers  with  zero  and  slightly 
favorable  pressure  gradients  has  been  made  where  data  obtained  at  Mach  numbers  up 
to  about  20  have  been  investigated.  The  data  were  classified  according  to  three  basic 
types  of  test  configurations;  namely:  (u)  flat  plate,  cone  and  hollow  cylinders;  (b)  two- 
dimensional  nozzle  wall;  and  (c)  axisjmimetric  nozzle  wall.  The  results  show  that  n 
is  primarily  a function  of  Reynolds  number  and  wall  temperature.  It  was  also  observed 
that  when  measurements  were  made  in  the  region  of  the  beginning  of  turbulent  flow 
(i.  e.,  the  end  of  transition  defined  as  the  point  where  occurs)  there  is  generally 
an  overshoot  in  the  parameter  n ranging  to  values  as  high  as  16.  As  a result  of  this 
survey,  it  appears  that  a relationshii)  exists  between  n and  Re,y  as  shown  in  Figure  10, 
where  it  is  noted  that  within  llie  overshoot  region  non-similarly  effects  are  prevalent. 

For  zero  mass  injection,  an  examination,  of  the  data  of  Reference  20  indicated  that 
the  axial  distribution  is  characterized  by  a value  of  n 9 at  the  forward  (transitional 
flow)  station  and  increases  to  some  unknown  value  (n  12)  at  the  end  of  transition.  A 
decrease  is  then  noted  to  values  of  the  order  of  7 at  points  sufficiently  far  removed 
from  the  end  of  transition.  This  behavior  is  characteristic  of  the  results  obtained  by 
Johnson -Bus hnelr'\  Moreovei*,  the  region  where  n decreases  with  increasing  Re^ 
persists  for  some  50  boundary  layer  thicknesses  downstream  of  the  point  where  q is  a 
local  maximum.  This  fact  is  significant  in  the  design  of  ground  tost  experiments  where 
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turbulent  separation  data  are  desired.  The  boundary  layer  profile  shape  can  have  a 
gross  impact  on  local  flow  characteristics  (such  as  separation  lengths)  relative  to  the 
separated  region.  As  a consequence  of  ground  test  ^cility  Reynolds  number  limita- 
tions, data  often  are  acquired  where  the  end  of  transition  is  sufficiently  close  to  the 
separated  flow  region.  These  data  must  be  viewed  as  optimistic,  in  that  separation 
lengths  could  be  minimized  due  to  the  large  concentration  of  high  energy  flow  near  the 
wall. 

The  effect  of  mass  transfer  on  the  turbulent  boimdary  layer  velocity  profile  tends 
to  reduce  its  fullness  (i.  e. , the  exponent  n is  reduced).  It  is  interesting  to  note  that 
the  axial  distribution  of  n with  blowing  retains  the  same  general  characteristics  as 
the  non-blowing  value  which  is  noted  by  Reference  20  shown  in  Figure  35. 

In  assessing  the  analytical  and  experimental  studies  on  viscous  layer  thickness  for 
turbulent  hypersonic  boundary  layers,  the  design  engineer  has  had  to  rely  on  semi- 
empirical  development  based  on  experimental  data.  This  is  a consequence  of  the  diffi- 
culty of  modeling  the  turbulent  shear  structure  in  three-dimensional  analyses  as  well 
as  time  consumption  which  is  inherent  in  codes  of  this  nature.  On  the  other  hand, 
recent  data  has  shown  that  engineering  codes  (integral  type)  have  a tendency  to  over- 
predict mass  transfer  effects  (with  q - O’ ).  Moreover,  the  inclusion  of  the  angle  of 
attack  effects  on  the  viscous  thicknesses  has  resorted  to  curve  fitting  of  limited  data 
while  virtually  no  information  has  been  available  for  the  coupling  of  mass  transfer  and 
angle  of  attack. 

The  technical  objective  of  the  present  study  is  to  examine  the  effects  of  mass  trans- 
fer and  angle  of  attack  on  the  viscous  turbulent  boundary  layer  parameters  6 , K*,  o , 
and  n.  The  data  obtained  from  the  SAMSO  STRKKT-G  study  (References  1 and  2)  will 
be  used  to  define  the  modifications  for  including  angle  of  attack  effects  into  existing 
engineering  codes  which  are  expressed  as  fonctions  of  mass  transfer  only.  The  anal- 
yses will  also  include  molecular  weight  effects. 
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4. 2 DATA  PRESENTATION 


4.  2. 1 Viscous  Layer  Thicknesses:  a - X = 0 


Inasmuch  as  the  empirical  correlations  developed  in  this  study  require  a standard 
for  normalization  (i.  e. , non-blowing,  zero  angle  of  attack  conditions),  the  current 
state  of  the  art  engineering  prediction  capability  will  be  assessed.  Essentially,  the 
development  of  this  work  is  classic  in  nature  and  employs  the  usual  boundary  layer 
assumptions.  An  equation  of  state  is  assumed  for  a constant  pressure  boundary  layer 
that  relates  density  to  temperature  (enthalpy).  The  Crocco  linear  velocity  tempera- 
ture expression  is  then  used  together  with  a power  law  velocity  distribution.  Bloom 
and  Martellucci^^  adopted  such  an  approach  where  a Dorodnitzyn  type  compressibility 
transformation  was  used  on  the  coordinates  with  the  velocity  power  law.  The  viscous 
layer  thicknesses  can  then  be  expressed  as 

r -1 

rl  _1 

*'0 

and 

Noting  the  equation  of  state  ( i>/  />g)  h^/h  and  the  velocity  power  law  u/ug  - ’fr\,  for 
n the  velocity  power  exponent,  one  can  write 
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The  linear  Crocco  relation  is  then  used  for  the  static  enthalpy,  giving 
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n ^ 1 ^ n 

(n  + 1)  (n  »-  2)  ^ (n  + 1)  hg  ^ (n  + 1)  (n  -f  2) 

(11)  together  with  the  velocity  power  law  reduces  to 

0 A = n/  [ (n  + 1)  (n  + 2)  ] 
o o ^ •* 

and  it  is  easily  shown  that 
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Equation 
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The  above  was  developed  in  Reference  32  for  A = B = unity.  However  Walker  and 
28 

Schumann  noted  a disparity  when  comparing  Eqs,  (17)  and  (18)  to  data  for  A = B = 
unity.  The  parameters  A and  B were  then  arbitrarily  placed  into  Eqs,  (17)  and  (18) 
with  the  kinetic  energy  term.  Walker  then  compared  equation  (17)  with  the  parameter 

OQ 

A to  the  numerical  results  of  Persh  , who  used  a power  law  velocity  distribution  to- 
gether with  the  Crocco  linear  velocity-temperature  relations,  from  the  integration  of 
the  expression 
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Equation  (17)  was  found  to  closely  approximate  the  numerical  results  for  A -•  1.29. 
However,  an  examination  of  available  data“  indicated  that  a value  of  A = 1. 69  showed 
aopreement  between  theory  and  data  to  within  10'', . 


If  one  considers  the  reduced  form  of  Eq.  (17)  from  the  Bloom-Martellucci  result, 
one  can  write 
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which  would  imply  that  the  Walker-Schumann  modification  can  be  expressed  as 

h /h  - 1 + A r 
aw  e 2 ® 

However,  from  a conservation  of  energy,  the  product  A r must  be  < one  which  is  not 

the  case  when  matching  the  data.  This  dilemma  cannot  be  explained  by  the  authors. 

An  inspection  of  Figures  36  and  37  shows,  for  n - 7 at  various  values  of  h /h  , 
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that  Eqs.  (17)  and  (18)  are  in  good  agreement  with  data  for  A = 1. 69  and  B = 1. 16 
(also  suggested  in  Reference  28)  including  the  recen  tly  acquired  data  of  References  1 
and  2. 

An  alternative  to  the  approach  used  by  Walker  and  Schumann  is  to  employ  the  mod- 
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ified  Crocco  relation  suggested  by  Kutateladze  and  Leont'ev  . Here,  the  authors 
choose  to  include  a non-similarity  parameter,  k . which  is  defined  as 


f (pressure  gradient,  mass  transfer,  enthalpy  potential) 


More  specifically,  J is  defined  by  Eqs.  (4)  and  (5)  of  Section  3.  If  o:.j  expresses  Eq. 
(7)  in  terms  of  the  static  enthalpy,  there  results 
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(19) 


Equation  (19)  can  be  substituted  into  Eqs.  (12)  and  (13)  and  the  normalized  viscous 
layer  thicknesses  become  {whf‘»'c  i is  assumed  as  a constant) 
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where  it  is  noted  that  for  (,  = unity,  Eqs.  (20)  and  (21)  reduce  to  Eqs,  (17)  and  (18). 


The  important  thing  to  note  here  is  that  the  parameter  ^ is  correctly  connected 

with  the  enthalpy  potential  term  (1.  e. , h /h  ) and  not  the  kinetic  energy  term.  More- 

w 

over,  an  examination  of  the  STREET-G  data  indicates  that  the  parameter  ^ can  hi 
used  to  model  the  viscous  layer  thicknesses.  Inasmuch  as  ? = ? (yA  ),  farther  ex- 
perimental information  will  be  required  to  actually  establish  i = ? ( n ) in  Eqs.  (12) 
and  (13).  However,  if  one  considers  a value  of  ^ and  n in  the  energy  overshoot 
region  (see,  for  example.  Figure  30)  good  agreement  is  noted.  Until  farther  informa- 
tion can  be  obtained  relative  to  the  non -similarity  parameter,  it  is  suggested  that  the 
Walker-Schumann  method  be  used  for  engineering  purposes. 

As  noted  in  the  preceding  development,  an  accurate  value  of  the  momentum  thick- 
ness will  be  required  in  order  to  determine  the  corresponding  viscous  thicknesses  Hq 
* 

and  f>Q  . Moreover  this  value  will  also  be  required  for  the  correlation  of  viscous 
layer  thicknesses  involving  mass  transfer  (which  will  follow).  Traditionally,  the 
momentum  thickness  has  been  determined  from  the  incompressible  momentum  equa- 
tion and  modified  to  include  compressible  effects  through  some  type  of  transformation. 
Moreover,  the  skin-friction  coefficient  is  replaced  by  a Blasius  type  solution  of  the 
coordinate  Reynolds  number  (which  implies  a power  law  velocity  distribution)  and  is 
subsequently  modified  to  include  compressibility  with  the  Eckert  reference  enthalpy 
method.  This  can  be  expressed  as 


classic  Blasius  equation^^,  which  will  be  generalized  to  include  a variable  power  law^  ? 


there  results 
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Here,  the  parameter  C is  a function  of  n also  (see  Reference  12  or  34).  The  boundary 
layer  thickness  can  be  determined  from  a conservation  of  momentum  which  allows  for 
a solution  to  the  skin  friction  coefficient.  The  flat  plate  momentum  equation  becomes 
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Integration  of  L'q,  (23)  yields 
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and  for  a sharp  cone,  the  above  becomes 
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where  the  Mangier  Factor  (MF)  is  given  by 

(l^n)  /(3+n) 
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It  is  important  to  note  that  the  above  formulation  did  not  allow  for  n = n(s)  which  has 
been  experimentally  observed  (Figure  10).  If  one  considers  the  classic  1/7  power  law, 
C has  the  value  8. 74  (which  has  since  been  modified  to  8. 56  to  agree  with  measured 
data)  and  the  parameter  K (n)  is  equal  to  0. 0296.  Table  V gives  the  values  of  C,  K, 
and  MF  (as  a function  of  n)  that  should  be  used  in  calculating  viscous  layer  thicknesses. 

4.2.2  Velocity  Power  Law  Exponent 

Data  representing  the  velocity  exponent  as  a function  of  Re^  are  shown  in  Figure 
10  at  two  free  stream  Reynolds  number  conditions.  Here,  one  notes  that  the  exponent 
n increases  through  the  transitional  boundary  layer  region  into  the  overshoot  region  as 
shown  in  References  16  and  17.  It  is  interesting  to  note  that  the  classic  1/7  distribu- 
tion (n  = 7)  is  never  reached  which  is  somewhat  surprising  when  considering  a 5-foot 
model.  This  leaves  one  to  speculate  on  the  validity  of  published  data  on  models  of  the 
order  of  one  foot  in  length  relative  to  fully-developed  turbulent  flow  characteristics. 

Figure  38  shows  the  attenuation  in  the  velocity  exponent  as  a function  of  mass 
transfer  and  molecular  weight.  The  data  of  STREET-G  are  shown  together  with  data 
from  References  19  and  35,  The  data  have  been  correlated  as  a function  of  the  blowing 
parameter  in  the  form 

4.2.3  Viscous  Layer  Thicloiesses  with  Mass  Transfer:  a = 0 

Figures  39  through  41  show  the  viscous  layer  thicknesses,  normalized  by  non- 
injection values,  as  a function  of  the  skin-friction  blowing  parameter.  The  data  rep- 
resent several  injection  rates  with  three  injectant  gases.  In  all  cases,  correlations 
were  developed  based  on  data  fitting.  It  is  noted  that  boundary  layer  thickness  growth 
displays  a slight  curvature  with  blowing  while  the  displacement  and  momentum  thick- 
ness growth  appears  to  be  linear.  Considering  the  various  molecular  weights  and 
injectant  rates,  the  data  show-  a well-ordered  behavior.  The  correlations  developed 
for  the  viscous  thicknesses  with  blowing  are 
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0.264  I 
1 + 0.03b, 


4-=  1 


-=1.0.41  (32) 

Figure  42  shows  the  viscous  layer  thicknesses  data  of  References  1 and  2 com- 
pared to  the  prediction  technique  developed  within.  The  method  of  prediction  consid- 
ers the  analytical  expressions  given  by  Eqs.  (17)  and  (18)  together  with  Eq.  (27)  and 
the  empirical  relations  represented  by  Eqs.  (30)  through  (32).  For  example,  the  dis- 
placement thickness  in  the  presence  of  blowing  can  be  written  as 
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Eq.  (31) 


It  should  be  noted  that  the  curves  shown  in  Figure  42  were  based  on  a 1/7  power  law 
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with  the  Walker-Schumann“  modification  to  the  Bloom-Martellucci  results.  Excel- 
lent agreement  is  noted  which  further  substantiates  the  choice  of  the  non-blowing  vis- 
cous layer  parameters. 

4.  2.  4 Viscous  Layer  Thicknesses  with  Angle  of  Attack  and  Blowing 


Several  attempts  were  made  to  correlate  the  viscous  layer  thicknesses  at  angle 
of  attack  conditions.  For  example,  the  boundary  layer  thickness  normalized  by  its 
zero  angle  of  attack  value  as  a function  of  a / 0 ^ (i.e.,  ratio  of  angle  of  attack  to  conical 
half-angle)  were  not  successful.  Moreover,  the  suggested  correlation  of  Copper  and 
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Shaw  (shown  in  Figure  43)  which  defines 
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was  also  attempted  with  no  further  improvement. 

Other  attempts  were  made  to  correlate  angle  of  attack  data.  Shown  in  Figures  44 
through  46  is  the  correlation  developed  in  this  study.  Here,  the  viscous  layer  thick- 
nesses are  shown  normalized  by  their  zero  angle  of  attack  value  as  a function  ® 
and  bj^.  Moreover,  the  data  shown  represent  various  injectant  gases.  Unlike  the  cor- 
relation attempts  with  a / as  the  independent  variable,  the  data  shown  in  Figures  44 
through  46  have  a tendency  to  coalesce  in  a more  orderly  fashion.  Inasmuch  as  the 
data  appears  to  be  characterized  by  a hyperbola,  a correlation  of  this  type  was  made 
and  is  shown  by  the  curves  in  the  figures.  In  each  case  the  curve  drawn  through  the 

data  is  the  hyperbola  for  which  the  ordinate  equals  the  reciprocal  of  the  quantity 
1/5 

(1  + a/e„  t b,,  /6),  that  is 

c M,  Q 

6(q)  (a)  i)  (o ) 

6(Q^0)’  6Mq-0)  ’ 0(0=0) 

With  the  exception  of  the  leeward  angle  of  attack  data,  satisfactory  agreement  is 
accorded. 
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4.3  DISCUSSION  OF  RElSiU.  TS 

In  the  work  statement  of  the  proposal  (GE-RESD  N72647  Rl)  concerning  profile 
shapes,  the  velocity  ex|)onent  was  to  be  examined  relative  to  the  effects  of  Reynolds 
number,  transition  location,  blowing,  and  angle  of  attack.  After  careful  deliberation, 
it  was  found  that  the  angle  of  attack  dependence  on  n did  not  show  any  trends.  Since 
no  pattern  could  be  ascertained,  no  further  work  was  considered.  Moreover,  inasmuch 
as  three  reporting  profile  stations  (in  general)  were  recorded,  no  definitive  conclusions 
could  be  made  relative  to  transition  location  effects  on  the  exponent  n. 

As  noted  in  the  development  of  the  zero  blowing  viscous  layer  thicknesses,  a 
dilemma  existed  concerning  the  arbitrary  choice  of  the  con.stants  A and  R of  Walker 
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and  Schumann  (see  Eqs.  (17)  and  (18)).  While  the  choice  of  the  value  A 
appears  to  violate  the  conservation  of  energy,  the  results  were  shown  to  be  in  good 
agreement  with  experimental  data.  On  the  other  hand,  one  can  arbitrarily  recast 
Eq.  (17)  to  read 

where  the  arbitrary  constant  is  now  a multiplicative  product  to  the  enthalpy  poten- 
tial term  (i.  e. , h /h  ).  One  will  note  that  the  above  is  much  in  the  same  spirit  as 
that  proposed  in  Reference  12.  When  Eq,  (34)  is  matched  with  Eq.  (17),  for  A = 1.69, 

n = 7,  and  h /h  = 1,  a value  of  the  parameter  \l^  of  1, 62  yields  results  to  within  3% 
w aw 

of  Eq.  (17)  for  Mg  > 5.  For  values  of  Mg  < 5,  departures  up  to  25%  are  observed  to 
values  of  Mg  = unity. 

4.4  CONCLUSIONS 

Turbulent  boundary  layer  characteristic  parameters  s,  f,*,  0 , and  n have  been 
examined  relative  to  the  influence  of  Reynolds  number,  angle  of  attack,  mass  transfer, 
and  molecular  weight  variations  of  the  injectant  gas.  It  was  found  that  the  state  of 
boundary  layer  development  was  sensitive  to  the  velocity  power  law  exponent  as  well 
as  the  viscous  layer  thicknesses.  Several  of  the  pertinent  observations  of  this  study 
are  as  follows: 

1.  The  velocity  power  law  exponent  showed  significant  departures  from  the 
classic  n = 7 value,  characteristic  of  fully-developed  flows,  to  values  of 

n >12  in  the  overshoot  region.  The  latter  defined  as  the  region  just  after 
transition.  As  a result,  the  reference  momentum  thickness  relations  had 
to  be  modified  to  include  a general  velocity  exponent. 

2.  The  prediction  capability  for  the  viscous  layer  thickness  parameters  has 
been  shown  to  yield  excellent  agreement  with  data  for  the  reference  case 

( ft  = X = 0)  where  a modified  Crocco  relation  was  used  with  a generalized 
velocity  power  law.  Specifically,  the  recommended  correlations  for  g, 
f>Q,  and  f>*  are  given  by  Eqs.  (27),  (17),  and  (18),  respectively,  where 
n = 7,  A = 1.  69,  and.B  = 1. 16. 
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3.  With  mass  transfer  effects,  the  velocity  power  law  exponent  attenuated  to 
values  of  n < 3.  On  the  other  hand,  the  viscous  layer  thicknesses  were 
augmented  above  the  non-injectant  values  as  might  be  expected.  Concerning 
the  latter,  correlations  were  developed  which  include  both  the  mass  trans- 
fer parameter  and  molecular  weight  ratio.  Prediction  capability  was  then 
afforded  by  combining  zero  injectira  analysis  with  the  mass  transfer  data 
correlations.  Specifically,  the  recommended  mass  transfer  correlations 
are  given  by  Eqs.  (30),  (31),  and  (32). 

4.  When  angle  of  attack  was  combined  with  blowing,  the  viscous  layer  thick- 
nesses did  not  appear  to  correlate  using  methods  previously  reported  in 
the  literature.  However,  a hyperbolic  function  was  foxmd  to  reasonably 
correlate  the  thicknesses  (when  normalized  by  their  a = 0 values)  as  a func- 
tion of  Of,  and  the  blowing  parameter  bj^  a . The  recommended  correlation 
is  expressed  by  Eq.  (33). 
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SECTION  5 

TURBULENT  HEAT  BLOCKAGE  AND  SKIN- FRICTION 
REDUCTION  DUE  TO  MASS  INJECTION 


5.1  BACKGROUND 

The  alleviation  of  the  high  heating  rates  encountered  by  surfaces  of  hypersonic 
vehicles  has  been  recognized  as  an  important  problem.  One  of  the  cooling  methods 
that  has  shown  oromise  is  mass  transfer  cooling,  wherein  a foreign  material  is  trans- 
ferred from  the  vehicle  surface  into  the  boundary  layer.  This  has  a two-fold  advantage 
in  attenuating  the  heat  transfer  problem.  The  transferred  coolant  may  absorb  heat  from 
the  boundary  layer  through  a phase  change  (sublimation,  evaporation,  melting,  etc.) 
as  well  as  providing  high  thermal  heat  capacitance.  In  addition,  it  has  been  shown  that 
the  introduction  of  a material  (with  its  normal  velocity  component)  at  the  surface  acts 
to  decelerate  the  flow  and,  consequently,  to  reduce  the  skin  friction.  This  also  implies 
a reduction  in  heat  transfer  at  the  wall. 

When  examining  the  available  literature  concerning  mass  transfer  experiments,  it 
is  noted  that  (1)  the  bulk  of  the  data  was  obtained  under  supersonic  flow  conditions 
(M  < 5)  and  (2)  the  model  geometries  consisted  primarily  of  flat  plates,  sharp  cones, 
and  cylinders  in  cross-flow,  precluding  pressure  gradient  effects.  Moreover,  a care- 
ful inspection  of  engineering  heat  transfer  prediction  techniques  indicates  that  empiri- 
cally derived  correlations,  developed  from  the  transpiration  experiments,  have  been 
employed  to  determine  the  level  of  heat  attenuation  resulting  from  blockage  by  the 
injectant  gas.  However,  examination  of  the  available  hypersonic  turbulent  boundary 
layer  data,  where  mass  transfer  effects  are  present  indicate  that  the  existing  empirical 
correlations^®*  ^®*  for  turbulent  boundary  layers  do  not  sufficiently  characterize  the 
heat  blockage  because  the  effects  of  Mach  number  and  wall  temperature  are  not  taken 
into  account  (see  Figure  47) . 

In  general  most  analytical  investigations  have  centered  around  the  flat-plate  or 
sharp  cone  geometry  to  allow  a more  tractable  mathematical  solution.  The  most 
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successful  approaches  involve  the  solution  of  the  momentum  integral  equation  coupled 
with  suitable  compressibility  transformations  and  Colburn's  analogy  to  predict  skin- 
friction  and  heat  transfer.  These  approaches  will  be  examined,  modified,  and  e}q)anded: 
first,  to  include  the  recently  acquired  experimental  data  of  Martellucci  et  al  ^ at 
hypersonic  (M  > 5)  conditions,  and  second,  to  account  for  the  effects  of  both  Mach 
number  and  wall  temperature,  with  the  objective  of  providing  improved  engineering 
design  relationships  suitable  for  inclusion  in  existing  semi-empirical  turbulent  boundary 
layer  codes. 

(a)  Heat  Blockage 

There  have  been  many  experiments  performed  on  various  geometric  shapes  in  low 
and  high  speed  tunnels  to  determine  the  effect  of  heat  blockage  on  a transpiration-cooled 
surface.  Experiments  of  this  nature  have  been  performed  to  simulate  ablation  effects 
experienced  on  re-entry  vehicles.  However,  analytical  predictions  have  weighed  heavily 
on  the  Reynolds  analogy  to  determine  the  attenuation  of  heat  flux  as  a consequence  of 
mass  injection,  thereby  requiring  recourse  to  experimental  data  and  subsequent  engi- 
neering correlations.  Of  the  three  empirical  correlations  alluded  to  previously 
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(Figure  47),  only  that  of  Walker**^  was  a direct  curve-fit  of  the  available  experimental 
data  at  that  time.  Lees'  expression  was  derived  analytically  based  on  a Couette-flow 
assumption.  Bartle  and  Leadon^^  correlated  their  e^qperimental  data  in  terms  of  an 
effectiveness  parameter.  They  suggested  that  the  heat  blockage  could  be  expressed  in 
terms  of  the  effectiveness  as  shown  in  Figure  47.  A comparison  of  these  correlations 
with  the  data  of  Martellucci  et  al  ^ (STREET-G)  at  M = 8 and  Pappas  & Okuno‘^2  ^t 
M ^ 0 shows  a definite  deficiency  in  accounting  for  the  effect  of  Mach  number. 

Two  of  the  most  pertinent  analytical  investigations  of  turbulent  boundary  layer  heat 
and  mass  transfer  employing  compressibility  transformations  are  those  of  Spalding, 
Auslander,  and  Sundaram'^3  and  Kutateladze  and  Leont'ev. ^2  Spalding  et  al  introduced 
compressibility  transformations  for  skin-friction,  heat  transfer  (Colburn  analogy)  and 
Reynolds  number  (viscous  effect) . The  developed  technique  utilized  all  the  available 
experimental  data  at  the  time  (zero  pressure  gradient,  supersonic  flow,  M < 5) . 
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Spalding  et  al  did  not  include  a Mach  number  or  wall  temperature  effect  in  their 
viscous  transformation  but  recognized  its  potential.  On  the  other  hand,  Kutateladze 
and  Leont'ev  included  the  effect  of  Mach  number  in  their  theory  but  neglected  the 
viscous  effect. 

Most  of  the  turbulent  experimental  heat  blockage  data  reported  up  to  the  present 
time  were  obtained  for  supersonic  flow  over  flat  plates  and  sharp  cones  (having 
negligible  pressure  gradients) . Bartle  and  Leadon'^^  and  Leadon  and  Scott^  reported 
data  for  heat  blockage  due  to  mass  injection  (air  and  foreign  gas)  for  a Mach  number 
range  of  Mg  = 2.0  to  3.2  with  corresponding  wall  temperature  ratios  (see  Table  VI). 
Pappas  and  Okuiio‘^2  reported  heat  blockage  data  with  air,  helium,  and  Freon  injection 
over  a sharp  cone  for  a Mach  number  range  of  Mg  = 0.7  to  4.35  and  =»  Ty.  Tewfik 
et  al  reported  subsonic  air  injection  data  over  a cylinder  in  cross-flow.  More 
recently  hypersonic  air  i.njection  data  on  a flat  plate  were  obtained  by  Dariberg46  at 
Mg  =6.7  which  indicated  a V'all  temperature  effect  on  heat  blockage. 

The  recent  STREET-G  data  of  Martellucci  et  al  ^ obtained  using  a sharp  porous 
cone  at  Me  = 7.1  (at  two  different  wall  temperatures  for  air,  helium,  argon,  and  Freon- 
12  injection)  supplies  the  bulk  of  the  hypersonic  flow  data  required  to  take  into  account 
the  effects  of  Mach  number,  wall  temperature  and  injectant  molecular  weight  on  heat 
blockage.  Extensive  angle  of  attack  heat  transfer  data  are  also  reported  in  References 
1 and  2. 

An  attempt  to  account  for  pressure  gradient  effects  was  made  by  Laganelli  et  al 
in  a porous  blunt  body  investigation;  however,  the  authors  noted  that  the  data  obtained 
were  questionable  due  to  the  sensitivity  of  the  tunnel  free-stream  conditions  and  surface 
sensor  reponse. 

(b)  Skin-Friction 

Relative  to  the  abundant  heat  blockage  literature,  there  is  a paucity  of  applicable 
skin-friction  experiments.  Only  Dershin  et  al  in  their  flat  plate  air-injection 
investigation  have  measured  local  skin-friction  reduction  in  a supersonic  flow  (Me=3.18) 
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environment.  Dershin's  data  tended  to  verify  Reynolds  analogy  when  correlated  with 
the  heat  blockage  data  of  Bartle  and  Leadon.^^  The  subsonic  flow  measurements  of 
local  skin  friction  reduction  with  air  injection  performed  by  Kendall  et  al  and 
Goodwin®®  serve  to  complement  the  supersonic  data  of  Dershin. 

Pappas  and  Okuno^l  reported  total  (averaged  over  the  body)  skin-friction  meas- 
urements over  a porous  sharp  ccne  for  the  Mach  number  range  of  0 to  4.3.  Data  was 
obtained  for  air,  helium  and  Freon-12  injection.  This  data  together  with  that  of 
Tendeland  & Okuno®^  (Mg  = 2.55)  exhibited  a significant  Mach  number  dependence 
causing  a disparity  with  the  local  measurements  of  Dershin  et  al.  Dershin  suggested 
that  the  difference  may  be  due  to  small  local  pressure  gradient  effects  that  become 
apparent  for  total  skin-friction  on  a cone  but  may  not  be  relevant  for  local  flat  plate 
skin-friction.  Walker^®  directly  correlated  the  total  (or  average)  supersonic  skin- 
friction  reduction  data  reviewed  above  as  a function  of  Mach  number  and  injectant 
molecular  weight  without  employing  a theoretical  base  as  did  Spalding  et  al.'^® 

In  assessing  the  experimental  investigations  of  heat  blockage  and  skin-friction 
reduction  in  turbulent  boundary  layers,  one  notes  that  the  bulk  of  the  previously 
reported  experimental  data  were  obtained  at  supersonic  (M  < 5)  flow  conditions.  Only 
the  recent  measurements  of  Danberg^®  and  Maitellucci  et  al  (STREET-G)  supply 
the  much  needed  hypersonic  (M  ' 5)  flow  information  required  to  critically  examine 
state-of-the-art  engineering  design  relationships.  In  addition  the  STREET-G  experi- 
ments provide  a wealth  of  heat  transfer  data  for  air  and  foreign  gas  injection  obtained 
at  angle  of  attack  to  be  used  in  validating  3-D  viscous  heating  design  codes.  The  data 
of  Pappas  & Okuno*^^'  ®1  and  others  together  with  the  STREET-G  data  supply  sufficient 
data  for  analysis  of  injectant  molecular  weight  influence  on  heat  blockage  and  skin- 
friction  reduction. 

A careful  examination  of  the  profile  data  obtained  in  the  STREET-G  study  indicates 
that  the  local  variations  across  the  boundary  layer  (e.g.,  velocity  and  temperature) 
have  not  reached  a fully-developed  characteristic.  Keeping  in  mind  that  the  model 
used  in  the  STREET-G  experiments  was  five  feet  long,  one  must  carefully  weigh  the 
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significance  of  "fully-developed  turbulent  characteristics".  Moreover,  the  credibility 
of  reported  turbulent  boundary  layer  data,  as  published  in  the  literature  on  much 
smaller  models  raises  many  questions.  In  addition,  the  assemblage  of  all  applicable 
data  illustrates  the  sensitivity  of  the  skin-friction  coefficient  and  Stanton  number  data 
to  the  character  of  the  flow.  If  the  flow  is  not  fully-developed,  the  question  arises  as 
to  what  is  the  proper  value  of  Sto  (or  Cf^)  for  normalization  purposes.  As  the  blowing 
rate  increases,  the  transitional  (overshoot)  region  moves  forward  along  the  surface  of 
the  cone  (or  plate) . The  problem  then  becomes  to  determine  whether  to  normalize  the 
blowing  data  with  Sto  (°^  same  wetted  length  or  other  distance  based  on 

some  equivalent  boundary  layer  thickness  parameter  (such  as  Re^ ). 

The  two  most  pertinent  analytical  investigations,  Spalding  et  al  and  Kutateladze 
and  Leont'ev^^  employ  compressibility  transformations  to  the  solution  of  the  momentum 
integral  equation.  The  viscous  transformation  employed  by  Spalding  et  al  was 
empirically  developed  from  the  supersonic  data  available  at  the  time.  The  result  v/as 
that  the  method  did  not  satisfactorily  take  into  account  the  effects  of  Mach  number  and 
wall  temperature  exhibited  by  the  inclusion  of  the  recently  acquired  hypersonic  data  of 
Danberg  and  Martellucci  et  al.  They  did,  however,  suggest  that  the  viscous  trans- 
formation could  be  modified  to  account  for  Me  and  T^y  as  more  data  becomes  available. 

Kutateladze  and  Leont'ev  did  not  include  a viscous  transformation  but  introduced 
the  concept  of  a critical  blowaway  parameter  based  on  an  infinite  Re^  - limiting  law. 
However,  their  derivation  resulted  in  critical  blowaway  values  decreasing  with  increas- 
ing Mach  number  which  is  contrary  to  the  observed  experimental  data. 

The  technical  objectives  of  the  present  study  are  to  examine,  modify  and  expand  the 
aforementioned  heat  blockage  and  skin-friction  reduction  correlations  (both  analytical 
and  experimental)  in  order  to  include  the  recently  acquired  data  of  Martellucci  et  al 
and  Danberg  obtained  at  hypersonic  flow  conditions  (M  - 5)  and  to  account  for  the 
effects  of  Mach  number,  wall  temperature  and  injectant  molecular  weight.  Ti.e 
resulting  semi-empirical  relations  are  then  to  be  employed  in  updating  the  GE-RESD 
3-D  Viscous  Code  to  predict  angle-of-attack  eifects  on  surface  heating.  In  addition. 
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this  study  serves  to  reduce  the  data  of  Martelluoci  et  al  to  a form  suitable  for  addi- 
tion to  the  general  literature  on  heat  blockage. 


5.2  ANALYSIS 


5.2.1  General  i^proach 


Considerable  success  in  correlating  heat  transfer  blockage  and  skin-friction 
reduction  as  a function  of  Mach  number  and  wall  temperature  was  obtained  by  adopting 
an  approach  similar  to  that  of  Spalding  et  al  which  employs  compressibility  trans- 
formations to  the  solution  of  the  momentum  integral  equation.  For  a cone  with  mass 
injection  one  may  write 


1 ± 
r ds 


(r  <?) 


Cf  Cf 

~ + X = - a+B^) 


(35) 


where  Bu  is  defined  as  2 A /Cf . 

For  a sharp  cone;  r = s sin  0 c»  and  the  momentiAm  integral  equation  in  terms  of 
Reynolds  number  becomes 


d(RegRep  ) 
dReg 


Cf 

R«s  T ^®u) 


(36) 


Spalding  et  al  introduced  the  following  transformations  in  their  method: 

(a)  Compressibility  transformation 

— : Fq  = (Isle,  T,/Te.  B^)  (37) 

Cf 

(b)  Viscous  transformation 


= Fo  = F^  (Mg.  'Iw/Tg.  Bu)  (38) 

= Fg  = Fg  (Me,  T^/Tg,  By)  (39) 

Reg 
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When  these  transformations  are  applied  to  the  momentum  integral  relation, 
Eq.  (36),  there  results 


d(Fg  Re^  . F^RCg) 


F Re 
s s 


2 


a +Bu) 


F F 
c s 


(40) 


Raiding's  basic  postulate  asserts  that  unique  relations  exist  between  (Fq  Cf/2), 
(F^  Re^)  and  (FgReg).  This  can  only  be  true  if,  in  Eq.  (40),  (1  + Bu)  Fg /(Fq Fg)  is 
equal  to  a constant  for  all  boundary  layers.  Consideration  of  the  pairticular  case  of  the 
uniform  density  zero  mass  transfer  boundary  layer  shows  that  this  constant  must  be 
equal  to  unity.  Therefore 


Fs  = ^ (1  +Bu) 


(41) 


Thus,  it  is  seen  that  Fg  is  not  indepeixient  of  the  Fg  and  F^  transformations. 
This  means  that  only  Fc  and  F^  need  be  deduced,  either  analytically  or  empirically, 
to  afford  a solution  of  Eq.  (40) . 

Rewriting  Eq.  (40)  yields 


d(Fg  Re<,  • Reg) 
d RCg 


^0  Cf 

= — Res  (Fc— ) (1  +Bg) 


(42) 


The  skin  friction  coefficient  for  an  incompressible  flat  plate  flow  can  be  expressed 
as 


or 


-m 


(43) 


Cf  j,  -m 

Fc  2 ~ 2 ^ 


where  E and  m are  constants  related  to  the  velocity  power  law  exponent;  specifically, 
E = 0.013  and  m =0.25  for  n = 7. 


Eq.  (42)  then  becomes 


d(FgRe0-  Re^) 
die 

s 


E F.  m+1 

7 F (Ffi  Ree  ‘ Keg)  Re^  a +B„) 
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(44) 


which  integrates  to  (assuming  constant) 


_ „ _ .m+1  m+2 

(Fg  Reg -Reg)  F.  Reg 

= _ _JL  n +B  \ — ~ 

m+1  2 F^  ' u'  m+2 

c 

Rearranging  yields 


(45) 


m 

m+1 


(46) 


Reintroducing  Eq.  (43)  yields  the  result 


El  - A 

2 ~ 2 [ 2 V m+2 


^ \ rt  . T3  X D 1 
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m 
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For  the  zero  blowing  case  (B^,  = 0),  Eq.  (47)  reduces  to 


^fo 


Co  2 


_E 
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L ( ilLtL^  ^ c 

2 \m+2/ 


m 
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which,  when  combined  with  Eq.  (47),  results  in  tiie  skin-friction  reduction  relation 


(49) 
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If  one  employs  the  use  of  Reynolds'  analogy  (St  = Cf/2)  at  this  point,  the  heat 
blockage  relation,  analogous  to  Eq.  (49),  becomes 


m m 

m+1  m+1 


(1  + Bjj) 


(50) 


where  is  defined  in  terms  of  the  Stanton  number: 

Bh  (61) 

5.2.2  Deduction  of  the  Transformation  Functions 

Spalding  et  al  adopted  a definition  for  Fg  which  has  been  commonly  accepted  in 
the  scientific  community,  namely; 


where  z is  the  velocity  ratio  u/ug. 

In  the  above,  an  equation  of  state  is  generally  invoked  together  with  the  boundary 
layer  assumption  of  zero  pressure  gradient  normal  to  the  surface.  The  linear  Crocco 
relation  is  then  applied  and  numerical  integration  is  performed  to  evaluate  Eq.  (52). 

In  order  to  simplify  the  complexity  of  the  above  integral,  it  will  be  assumed  that  Fg 
can  be  approximated  by  the  expression 


Fc 


‘ _ l!L  ("  dz 


(53) 


where  the  starred  quantity  is  based  v . che  Eckert  reference  enthalpy;  then 


(54) 
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Note  that  in  the  limit  (using  L’Hopital's  rule)  as  -*  0 


which  is  the  transformation  function  for  a non-injection  case^^. 

A comparison  of  the  Fq*  approximation  with  the  numerically  exact  solution 
from  the  tables  of  Spalding  et  al  is  given  in  Figure  48  for  a realistic  range  of  Mach 
numbers  (0  < Me  ^ 8)  and  wall  temperature  ratios  (1  < T^/Tg  < 8) . The  comparison 
indicates  excellent  agreement  In  the  supersonic  regime  (M<  5)  but  a significant  vari- 
ation with  increasing  values  of  Mach  number;  however,  this  is  to  be  expected  as  seen, 
for  example,  in  a turbulent  heating  fli^t  test  analysis  by  Laganelli^'^.  The  important 
thing  to  keep  in  mind  here  is  that  this  approximation  is  purely  analytical  while  the  F^ 
viscous  transformation  function  (as  will  be  seen)  is  deduced  from  e^qierimental  data. 
The  net  effect  of  both  of  these  transformations  will  satisfactorily  account  for  the 
hypersonic  situation. 

While  most  authors  agreed  on  the  analytical  form  for  the  F^,  transformation, 
significant  disagreement  prevailed  in  the  formulation  of  the  F^  function.  Since  there 
existed  considerable  uncertainty  about  the  nature  of  F^ , Spalding  et  al  deduced  the 
following  empirical  relation  (based  on  all  the  air  injection  data  available  at  the  time 
prior  to  1962) 


F,  =il^(l+B^)  (56) 

n-w 

which  does  not  contain  a Mach  number  effect. 

Consequent  modification  of  this  relation  to  include  the  experimental  air  injection 
data  of  Martellucci  et  al.  Danberg'^®  and  Pappas  & Okuno'^2  along  with  the  Eckert 
reference  state  concept  resulted  in  the  empirical  relation 


F„  = 


— a + Bu) 
<1  ^ 


(57) 
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where 

= (Tw/Te)  ® + I ^*e  (58) 

Eq.  (58)  takes  into  account  the  effects  of  Mach  number  and  wall  temperature  through 
a purely  empirical  correlation  of  ejqjerimental  data  and  judicious  use  of  the  Fq  tables 
of  Spalding  et  al.  Note,  again,  that  in  the  limit  as  Bu  0,  Eq.  (57)  reduces  to 


F 


0 


o 


(59) 


which  is  the  correct  zero  injection  viscous  transformation  using  the  Eckert  approach®^. 

Substitution  of  Eqs.  (54),  (55),  (57),  and  ^9)  into  the  skin  friction  e:iq)ression,  Eq. 
(49),  results  in  (for  m = 1/4) 


(^)  ■ i{^-)  ('•■.) 

KGq  -J 


u;-  1 

5 


(60) 


where  it  is  noted  that  in  the  limit  as  B,j  -♦  0,  Cf/C[^  approaches  unity  as  required. 
The  anajlogous  equation  for  Stanton  number  is  simply  (by  Reynolds  analogy) 


Bh 


n 8 


<*-•-1 


(61) 


\^here  Bj^  has  replaced  B^^. 

5.2.3  Verification  With  Exi^erimental  Data 


A comparison  of  predictions  using  Eq.  ((51)  and  a wide  x’ange  of  air  injection 
experimental  data  is  given  in  Figure  49  where  the  heat  blockage  (represented  by  the 
Stanton  number  ratio,  St/St^)  is  plotted  as  a function  of  the  blowing  parameter  x /Stg. 
All  of  the  STREET-G  ^ data  points  were  obtained  at  the  high  Reynolds  number  (Re„o  /ft. 
= 3. 7 X IqO)  . iTie  Danberg  data  points  at  'l\v/Tg  “ 'i*  1 are  not  the  exact  St/Sto  values 
reported  in  Reference  46.  Tiie  (iata  points  in  Figure  49  were  obtained  by  normalizing 
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his  reported  St  values  with  a second  set  of  StQ  values  he  reported  at  the  same  condi- 
tions. It  was  felt  that  the  second  set  which  consisted  of  a larger  number  of  data  points 
would  prove  more  reliable.  This  conclusion  seems  to  be  borne  out  by  the  correlation 
curve  (Me  = 7,  T^/Tg  = 4)  in  Figure  49.  The  Pappas  and  Okuno^^  data  plotted  in 
Figure  49  are  their  "thermocouple  no.  7"  data  which  were  measured  at  their  highest 
reported  Reynolds  numbers.  Further  commc.nts  concerning  the  choice  of  the  data 
included  above  will  be  reflected  upon  in  the  Results  Section. 

The  wall  temperature  effects  were  taken  into  account  (Eq.  (58))  through  the  use 

AO 

of  the  Fq  transformation  function  tables  of  Spalding  et  al.  As  seen  from  the  figure 
the  correlation  curves  bracket  the  data  very  well,  verifying  the  approximations  for 
Fq  (Eq.  (54))  and  F^  (Eq.  (57)). 

The  influence  of  injectant  molecular  wei^t  on  heat  blockage  has  been  taken  into 
account  utilizing  the  foreign  gas  injection  data  of  Pappas  and  Okuno'^^,  Leadon  and 
Scott44  a:id  Martellucci  et  al  1*  ^ (STREET-G) . Correlation  of  heat  blockage  with 
foreign  gas  injection  is  presented  in  Figure  50  where  data  for  helium,  argon,  and 
Freon-12  injection  are  compared  with  the  correlation  of  Eq.  (61)  with  the  blowing 
parameter  Bh  redefined  as 


A 

■gr 


(62) 


Considering  the  data  scatter,  Eq.  (61)  does  an  excellent  job  of  correlating  the  a\^il- 
able  data. 


When  the  abscissa  oi  Figure  50,  (a/SIq),  is  multiplied  by  the  specific  heat  ratio 
(Cpin./Cp  j ^ the  result  is  Figure  51.  Attempts  to  correlate  the  data  using  molecular 
wei^t  as  a parameter  were  unsuccessful.  However,  the  specific  heat  factor  did 
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collapse  the  data  to  a band  reasonably  predicted  by  the  correlation.  The  Cp  values 


employed  for  the  various  Injectants  were: 

Ihjectant 

Molecular  Weight 

Cp(BTU/lb  °F) 

Air 

28.966 

.2401 

Argon 

39.944 

.1244 

Freon-12 

120.93 

. 1458 

Helium 

4.003 

1.2413 

5,2.4  Zero-blowing  Predictions 

Consistent  wltli  this  development  is  the  use  of  the  Kckert  reference  enthalpy  con- 
cept in  accounting  for  compressibility  effects  on  the  zero-blowing  skin-friction  and 
Stanton  number  values  used  as  the  nornuilizatlon  terms  in  P'q.'s.  (60)  and  (61),  These 

expressions  (based  on  the  velocity  power-law  concept  for  flat  plate  fully-developed 
54 

turbulent  flows)  are 


1 


= ,0296  (MF)  f j Reg  ^ 

(63) 

/q  _ 1 /r. 

= .0296  (MF)  f Pr  Re..  ' 

(64) 

1/5 

where  MF  is  the  Mangier  factor,  (2.25)  - 1.  176  for  a sharp  cone,  and  the  compressi- 

bility term 


(65) 


is  totally  consistent  with  the  preceding  analysis  (see  Fq.*s  (4!)),  (55)  and  (60)  with  m - 
2/(n+l)  = 1/4  for  n --  7,  the  classic  fully -(level  oped  turbulent  velocity  power  law  expon- 
ent). The  STRFPM'-G  "cold  wall"  zero-blowing  heat  transfer  data  are  plotted  in  Figure 
52  as  a function  of  wetted-length  Reynolds  number.  Two  groups  of  local  Stanton  number 
data  are  included  to  Illustrate  data  repeatability.  The  data  are  compared  with  the  pre- 
diction of  Fq.  (61)  which  is  based  on  the  one-seventh  power-law  profile  approximation 
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and  Colburn's  analogy,  iispectlon  of  Figure  52  shows  the  transitional  overshoot 
region  prior  to  fully-developed  turbulent  flow.  Excellent  agreement  is  observed  be- 
tween the  prediction  and  the  data  for  the  fully-developed  region. 


5o3  RESULTS 


5.3.1  Mass  Transfer  Correlations 


llie  skin-friction  reduction  relation.  Eq.  (60).  and  the  heat'blockage  expression. 
Eq.  (61).  that  have  been  developed  are  expressed  in  the  implicit  parameters  and 
B^.  which  contain  the  terms  Cf  and  St,  not  Cf^  and  St^.  As  engineering  design  corre- 
lations, these  expressions  would  require  Iterative  solutions  In  order  to  obtain  the  de- 
sired result.  Therefore,  more  tractable  correlations  for  design  purposes  should  be 
expressed  in  terms  of  explicit  blowing  parameters  (l.e. . terms  containing  the  a priori 
known  values  of  Cj^or  Sto). 

This  was  accomplished  through  the  use  of  the  critical  blowaway  parameter,  bcr, 

12 

introduced  by  Kutateladze  and  Leont'ev.  They  defined  b^j,  as  the  limiting  value  of 
the  explicit  blowing  parameter  when  Cf/Cf^  or  St/Stp  vanishes,  T.>e  explicit  blowing 
parameters  with  provision  made  for  foreign  gas  injection  (analogous  to  Bu  and  Bh)  are 
defined  as 


2>. 

Cf 


bh 


(66) 


(67) 


The  critical  blowaway  parameter  bcr  was  found  to  be  related  to  Mach  number  and 
wall  temperature  through  the  u>  parameter  defined  in  Equation  (58).  The  expression 
which  relates  b^ j.  to  oj  is  given  as 

bgj.  = exp  r 1.676  (<*>+ 0. 161) 
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(68) 


This  result,  Eq.  (68),  represents  the  approximation  for  b^j.  that  provides  the 
most  suitable  agreement  with  the  correlation  Eqs.  (60)  and  (61)  developed  in  the 
Analysis  Section.  Equation  (68)  is  plotted  in  Figure  53  and  compared  with  all  the 
available  heat  transfer  and  skin  friction  air  injection  data  (See  Table  VI).  Each 
symbol  represents  the  average  value  of  b^j^.  extracted  from  each  data  source.  The 
massive  blowing  data  of  Fogaroli  and  Saydah*^^  was  included  for  comparison  at  higher 
values  of  bcj..  The  incompressible  value  for  b^j.  (i.e.  when  = 1)  was  determined 
to  be  bgj.  = 7 in  order  to  match  the  predictions  of  Eqs.  (60)  and  (61).  This  value  is 
at  variance  with  the  reported  incompressible  value  of  b^j,  ^ 4 as  derived  by  Kuta- 
teladze  and  Leont'ev.  However,  their  result  (from  the  definition  below) 


bcr 


1 + 2 (y/6  ) 
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(69) 


was  based  on  the  assumption  that 


1 2(y/ A ) >•  1 (70) 

If,  instead,  an  effective  or  weighted  value  for  (y/.s  ) was  employed  in  Eq.  (69),  a 
larger  value  of  b^j.  would  result.  For  ex.'imple,  if  an  effective  value  of  (y/b  ) = 1/3 
were  assumed,  a value  for  bd-  of  approximately  7 would  result. 

Utilizing  this  concept  of  bgj.»  the  skin -friction  reduction  and  heat  blockage  equations 
may  be  expressed  in  terms  of  the  explicit  blowing  parmneters  b^  and  bj,  as  follows: 


The  form  of  these  expressions  was  inspired  by  a correlation  of  a series  of  exact  anal3^i- 
cal  solutions  (for  a wide  range  of  Mach  nuinlier  and  blowing  rate)  performed  by 
Kutateladze  and  liCont'ev.^® 


Li  Figure  54  the  heat  blockage  correlation,  Eq.  (72),  is  compared  with  the  air 
injection  data  reviewed  previously  in  the  discussion  of  Figure  49.  The  curves  gen- 
erated from  Eq.  (72)  provide  slightly  better  agreement  with  the  data  than  do  the 
equivalent  curves  generated  from  Eq.  (61).  Overall,  excellent  agreement  is  ob- 
served for  the  air  injection  comparison  over  the  Mach  number  range  0 < Mg  < 7 
for  three  different  geometrical  surfaces. 

A comparison  of  the  skln-frictlon  reduction  correlation,  Eq.  (71),  with  all 
available  air  Injection  data  Is  presented  in  Figure  55.  The  data  of  Dershln  et  al'^®, 
Kendall  et  al^9,  and  Goodwin®®  represent  direct  measurements  of  local  flat  plate 
skin  friction.  The  data  of  Pappas  and  Okuno®^  and  Tendeland  and  Okuno®^  repre- 
sent total  (averaged  over  the  body)  skin  friction  measurements  on  a cone  which  are 
denoted  by  Cp  as  distinguished  from  Cf.  Note  the  disparity  between  the  local  skin- 
friction  data  of  Dershln  et  al  (Mg  = 3. 18)  and  the  total  or  average  skin -friction  data 
of  Pappas  and  Okuno  (Mg  = 3.21).  Dershln  suggests  that  the  difference  may  be  due 
to  small  local  pressure  gradient  effects  that  become  apparent  for  total  skln-frictlon 
on  a cone  but  may  not  be  relevant  for  local  flat  plate  skin  friction,  to  any  event  the 
curves  generated  from  Eq.  (71)  for  the  range  0 < Mg  < 5 are  In  reasonable  agree- 
ment with  the  reported  data,  considering  the  disparity  In  the  data  itself.  The  curves 
were  generated  assuming  since  the  data  sources  did  not  report  the  ex- 

perimental wall  temperatures. 

Foreign  gas  Injection  data  for  b^^  (molecular  weight  influence)  were  added  to 
the  air  Injection  data  depicted  In  Figure  53  with  Figure  56  as  the  result.  Data  for 
helium  and  Freon-12  Injection  over  a wide  range  of  Mach  numbers  were  included 
with  relatively  Insignificant  scatter  around  the  air  injection  data  and  the  correla- 
tion curve,  Eq.  (68). 

Finally,  the  correlations  for  heat  blockage  and  skin-friction  reduction  with 
foreign  gas  injection  are  illustrated  in  Figures  57  and  58  where  St/St^  and  Cf/Cfo 
are  plotted  as  a function  of  bj,/bgj.  and  by/bg^.,  respectively.  Most  of  the  data  for 
air,  helium  and  Freon- 12  Injection  are  bracketed  by  the  correlation  curves  for 
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= 1 (subsonic  flow)  and  u = 2 (hypersonic  flow)  which  represent  a realistic  flight 
condition  range.  Some  of  the  helium  injection  data  fall  outside  the  subsonic  limit 
(oj  = 1)  in  both  figures.  Overall,  however,  the  two  Figures  57  and  58  verify  the 
success  of  the  correlation  equations  derived  for  heat  blockage  (F.q.  (72)  ) and  skin 
friction  reduction  (Eq.  (71)  ). 

5.3.2  Angle  of  Attack  Predictions 


The  heat  blockage  and  skin  friction  reduction  correlations  developed  in  this  study 
have  taken  into  account  the  effects  of  Mach  number,  wall  temperature  and  injectant 
molecular  weight  (foreign  gas  injection).  'I’hese  correlations  have  been  incorporated 
into  the  GE-RESD  3-D  Viscous  Code^^  in  order  to  generate  viscous  flow  solutions  for 
angle  of  attack  comparisons  with  applicable  STREET-G^  data. 

The  Three-Dimensional  Viscous  Code  applies  an  integral  boundary  layer  type 
solution  for  heating  and  wjill  shear  stress  to  sharp  cone,  sphere-cone,  and  bent  sphere- 
cone,  configurations  at  angle  of  attack.  The  integral  boundary  layer  solutions  are 
carried  out  along  streamlines  calcul.ated  in  tlu?  inviscid  sense  from  an  imposed  pres- 
sure field,  iiio  small  cross -flow  approxin\atic»n  is  made,  and  the  effects  of  angle  of 
attack  on  heat  transfer  arise  from  the  streamline  divergence  effect.  'ITils  is  included 
in  the  calculation  as  a scale  factor  or  metric  whicli  describes  the  rate  of  spread,  and 
which  at  zero  angle  of  attack  reduces  to  tlic  local  body  radius, 

Tlie  effect  of  angle  of  attack  on  surface  heat  transfer  was  exiunined  by  comparing 
solutions  from  the  GE  3-D  \iscous  Code  witli  .S'lHEET-G  cold  wall  ('I\^,  = 580  °R)  tur- 
bulent he.ating  data  obtained  on  the  windward  meridian  of  a sharp  cone  (5“  half-angle) 
at  Me  = 7. 1.  'I'he  surface  lieat  transfer  rale  as  a function  of  wotted  lengtli  along  the 
cone  for  three  values  of  angle  of  attack  ( <>  - (i’’,5“,  and  10'*)  is  presented  in  Figure 
59  (zero  injection)  and  Figure  00  (injection  rale  - .0015).  A comparison  of  the 
code  prediction  with  the  STRI.E  l -G  data  .-show.';  excellent  agreement  for  the  cases  with 
and  without  injection.  The  slight  disagreement  in  the  zei’o  injection  case  at  a - 10® 
(Figure  59)  is  offset  by  the  reia.iHsable  ju’edictive  capability  at  the  same  angle  of  attack 
case  with  injection  (Figure  OO). 
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5.3.3  Further  Comments 


(a)  Reynolds  Analogy 

Since  the  STREET-G  program  provided  both  surface  heat  transfer  data  and 
boundary  layer  profile  measurements,  with  and  without  blowing,  an  attempt  was  made 
to  Investigate  the  validity  of  the  Reynolds  analogy  experimentally  for  hypersonic  flow. 
Here,  directly  measured  Stanton  number  data  were  to  be  compared  with  skin  friction 
data  derived  from  the  profile  measurements.  The  skin-friction  coefficient  was  to  be 
determined  from  the  momentum  Integral  equation  (for  a sharp  cone) 


(73) 


where  the  blowing  rate  and  momentum  tliickness  were  measured.  Since  the  momentum 
thickness  was  recorded  at  three  stations  along  the  cone  surface,  the  slope  (do  /ds)  was 
obtained  from  a plot  of  O as  a function  of  the  wetted  length,  s. 

However,  because  of  the  low  enthalpy  potential  provided  by  the  AEDC  facility,  all 
the  terms  on  the  light  hand  side  of  Eq.  (73)  were  found  to  be  of  the  same  order  of 
magnitude,  requiring  subtraction  of  numbers  of  the  order  of  10"^.  This  numerical 
sensitivity  rendered  the  acquisition  of  meaningful  skin  friction  data  Impossible.  As  a 
consequence,  the  attempt  to  experimentally  verify  Reynolds  analogy  at  hypersonic 
conditions  had  to  be  abandoned. 

(b)  Concept  of  Fully-Developed  Flow 

The  STREET-G  experiments  have  demonstrated  that  a fully -developed  turbulent 
boundary  layer  did  not  exist  (from  a volocify  profile  power  low  exponent  point  of  view; 
i.e.  for  n = 7)  on  a five-foot  model.  This  point  has  been  discussed  extensively  in 
Sections  3 and  4 of  this  report.  This  knowledge  induces  one  to  speculate  on  the  validity 
of  experiments  previously  performed  using  models  on  the  order  of  one  foot  in  length. 
Several  attempts  were  made  to  correlate  the  STREET-G  heat  transfer  data  in  terms 
of  Reynolds  numbers  based  on  eciuivalent  momentun;  thickness,  enthalpy  deficit  thick- 
ness and  wetted  length  along  the  surface.  It  was  found  that  the  last  choice  contributed 


the  best  behavior  in  the  correlation  for  boat  transfer  attenuation.  However,  only  tlie 
data  toward  the  end  of  the  model  were  umphasi/.ed  since  they  were  the  least  likely  to 


be  influenced  by  the  velocity  exponent  overshoot  region.  It  was  for  this  reason  that 
only  the  data  reported  at  the  last  station  (and  highest  Reynolds  number)  of  the  models 
of  other  experimental  investigators  were  Incliided  for  use  in  this  study. 

5.4  CONCLUSIONS 

The  effect  of  Mach  number,  wall  temperature  and  injectant  molecular  weight  on 
heat  blockage  and  skin  friction  reduction  resulting  from  mass  injection  has  been  evalU' 
ated  for  turbulent  flow  over  surfaces  witli  negligible  axljil  pressure  gradients.  The 
evaluation  resulted  in  semi  -empirical  correlations  which  are  based  on  a modification 
of  a compressibility  transformation  theory  expa.nded  to  include  the  latest  hypersonic 
data  up  to  a free  stream  Mach  number  of  s.o.  Angle  of  attack  effects  on  heat  trans- 
fer were  examined  utilizing  Uie  mass  trimsfer  correlations  de\’eloped  in  this  sbady. 

Several  of  the  nuna;  Importanl  conclusions  dcxdved  from  this  study  are  as  follows: 

1.  It  has  been  shown  that  heat  transfer  attenuation  due  to  mass  injection  is 
heavily  influenced  l)V  Mach  number  and,  to  a lesser  degree,  by  wall 
temperature. 

2.  .A  semi-empirical  technuiuc  uas  developed  Uial  predicts  heat  blockage 
and  skin  fx’iction  reduction  due  lo  mass  imeclion  as  a function  of  Mach 
number  .'ind  wall  tcmper:*ture.  I'liis  method  was  ixased  on  an  analytically 
derived  compressibility  ir.uisform.atiou  .uid  an  empirically  derived  (based 
on  experimental  data)  vIj'Cohs  tra.'slonuation.  A critical  blowaway 
parameter  concept  was  i'i:i|)i(.>vt.-d  lo  provide  a m<n'e  tractable  set  of 
correlations  for  engineering  ocsign. 

3.  The  an.'dyticallv  dovixed  v'.impr.;ssii)iiitv  tr.-insformation  made  use  of 
the  Lckert  refcr-Mico  ctdhalpv  concept  with  consicleral)lc  success,  while 
the  viscous  transformation  development  met  with  sueccss  aided  by  the 
suggestions  of  .Spalding  ei  ;.!  a:u!  the  inclusion  of  recently  acquired 
hypersonic  flow  d;ita. 

•1.  The  resultant  heat  blucUage  and  skin  friction  reduction  correlations  pro- 
duced in  this  study  wen*  vc*rified  couclusivelv  by  comparison  with  data 
over  a wide  range  ot  .\I.ich  num!ie,i  (d  ■ 7)  and  wall  temperature 

(1  • 1\^. ■>'/.  riic  I'ffeets  of  iiijeclant  molecular  weight  were  ade- 

quately .accounted  for  i>y  a .-.pecifiv-  he.it  ratio  (C’pjj^.  C ) factored  into 
the  lilowing  jjar.anuU.  r.  i tu-  • voinmeiidid  correlations  for  Cf  C'f  and 


St/Stg  are  expressed  as  Eqs.  (71)  and  (72)  where  the  blowing  parameters 
are  defined  by  Eqs.  (66)  and  (67),  respectively.  The  exponent  w and  the 
parameter  b^j.  are  determined  from  Eqs.  (58)  and  (68),  respectively. 

The  use  of  the  Eckert  reference  enthalpy  concept  and  Colburn's  analogy 
resulted  in  satisfactory  prediction  of  the  zero  injection  STREET-G  heat 
transfer  data. 

Angle  of  attack  heat  transfer  was  satisfactorily  predicted  (when  compared 
witli  STREET-G  data)  by  the  GE-RESD  3-D  Viscous  Code  which  utilized 
the  heat  blockage  and  skin  friction  reduction  correlations  developed  in 
this  study. 
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If  one  equates  Equations  (A-4)  and  (A-6)  and  expands  the  total  derivatives  dT  and 
dP,  there  results 
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Considering  the  enthalpy  dependence  on  the  coordinates,  one  has 
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However,  from  the  boundary  layer  assumptions,  dP/d  y - 0;  hence. 
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where  it  is  understood  that 
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If  Equation  (A-11)  is  satisfied.  Equations  (A-3)  and  (A-10)  yield 
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Thus,  t last  term  on  the  RHS  of  Equation  (A-2)  becomes 
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Multiplying  Equation  (A-1)  by  u and  rewriting  Equation  (A-2)  using  the  definition 

Pr  = /j  C /k  and  h = CpT  and  finally  adding  the  resulting  two  equations,  one  has 
P 
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Noting  the  relationship 


and  unit  Prandtl  number.  Kq.  (A- M)  reduces  to 


which  has  the  solution 
2 

h — = constant 

The  constant  can  be  evaluated  using  the  wall  condition  u - 0,  such  that 

h = constant 
w 

and  we  can  write 
0 

h - h 

2 

The  enthalpy  gradient  at  the  wall  !)eeon\es 


Using  Uq.  (A-10),  we  can  write 


which  implies  an  adialiatic  wall  condition.  Thus,  the  arbitrary  constant  becomes 


h ^ h , ;md  Kq.  (A-Ui)  Is  written  as 

»!•  OtI*  * 


It  should  be  noted  tliat  Eq.  (A- 17)  was  developed  without  restricting  the  pressure 
gradient  to  zero.  Moreover,  an  examiriation  of  Eq.  (A-9)  indicates  that  for  a non- 
zero pressure  gradient,  tlie  temperature  is  not  necessarily  constant  along  the  surface 
even  though  the  enthalpy  must  be  constant  along  the  ^^urface.  On  the  other  hand,  if  the 
temperatures  are  fairly  low,  the  function  (6h/dp),,,  will  be  essentially  zero  such  that 
from  Eq.  (A-9),  regardless  of  the  pressure  gradient,  the  surface  temperature  must 
also  be  constant.  Now  at  high  temperatures,  the  corresponding  surfa.ce  condition  does 
not  have  to  be  constant  for  non-zero  pressure  gradients  but  must  be  constant  for 
zero  pressure  gradient.  Finally,  it  should  be  kept  in  mind  tnat  Eq.  (A-17)  is  valid 
only  for  unit  Prandtl  number. 


For  generality,  Eq.  (A-1)  can  be  multiplied  by  a constant,  say,  c,  which  yields 
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Combining  Lqs.  (A-18)  and  (A- 15)  there  results 
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which,  for  zero  pressure  gradient  conditions,  reduces  to 
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As  before,  one  possible  solution  to  the  above  is 
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where  at  the  wall,  u - 0,  one  has 


C' 


and  (flh/()y)  = -c(()u/,!y)  ^ 0.  Employing  the  conditions  at  the  edge  of  the  viscous 

w w 

layer,  there  results 
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(A-21) 

Inasmuch  as  the  derivation  of  Eq.  (A -21)  made  use  of  Eq,  (A- 15),  the  same  con- 
ditions apply  to  Eq.  (A-21);  namely,  that  Pr  = unity.  Moreover,  if  Pr  = unity,  the 
recovery  factor,  defined  as 
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must  also  be  unity.  Hence 
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The  static  enthalpy  can  then  be  expressed  as 
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and,  noting  the  definition  h 

aw 

distribution  is  written  as 


for  unit  Pra.idtl  number,  the  static  enthalpy 
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In  terms  of  total  enthalpy,  the  above  reduces  to  the  familiar  classic  linear  Crocco 
temperature-velocity  relationship,  namely 

I - h , 
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It  should  be  noted  that  the  wall  enthalpy  must  be  constant,  and,  since  the  pressure 
gradient  is  zero,  it  is  seen  from  Flq.  (A -9)  that  the  wall  temperature  must  also  be 
constant.  Thus,  the  conditions  that  must  be  satisfied  for  Eq.  (A -24)  or  (A-25)  to  be 
valid  are: 

(i)  unit  Prandtl  number 

(ii)  constant  pressure  (zero  pressure  gradient) 

(iii)  constant  surface  temperature 

In  Eq.  (A-25)  one  notes  a linear  relationship  between  H and  u which  implies  the 
existence  of  a local  Reynolds  analogy.  Moreover,  the  boundary  conditions  used  to 
obtain  Eq.  (A-25)  also  imply  that  the  relationship  is  valid  across  the  entire  boundary 
layer.  It  is  also  noted  that  mass  transfer  is  not  an  explicit  function  in  Eqs.  (A-1)  and 
(A-2),  and  hence  Eq.  (A-25)  as  well,  and  would  enter  only  through  the  boundary  con- 
ditions. 


As  a consequence  of  the  kinetic  energy  contribution  in  high  speed  flow,  Eq,  (A-14) 
can  be  re-expressed  as 
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Using  the  definition  of  total  enthalpy.  H h ‘ , the  above  becomes 
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For  gases,  the  Prandtl  number  is  close  to  unity  which  is  a basic  condition  for  the 

existence  of  similarity  between  the  hydrodynamic  and  thermal  profiles  (i.e. , Eq. 

(A-25)).  However,  departures  from  similarity,  which  are  a consequence  of  the  kin- 

.,2 

etic  energy  term  (Pr-1)  , can  be  accounted  for  by  considering  the  work  of 

Kutateladze  and  Ixjont'ev.  In  their  work,  the  authors  suggested  a relationship  of  the 
form 


(A-27) 


where 
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The  parameters  (.  and  r (recovery  factor)  are  both  functions  of  the  coordinate  and 
depend,  in  general,  on  the  pressure  gradient,  the  temperature  difference  imposed, 
and  the  magnitude  of  mass  transfer.  Spalding  noted  that  one  superfluous  function  had 
been  introduced  in  the  definitions  of  Eqs.  (A-27)  and  (A -28).  If  one  defines  the  re- 
covery value  as  a constant,  then  Eq.  (A-27)  together  with  Eq.  (A-28),  defines  the 

parameter  ? ( c)*  If.  the  other  hand,  r is  allowed  to  vary  with  (, , then  the  defin- 

0 

itions  are  incomplete. 


By  considering  the  recovery  factor  to  be  constant,  Eq,  (A-27)  can  be  re-written 
in  terms  of  the  static  and  total  temperature  (Crocco  variables)  as 
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then  the  total  temperature  becomes  (Crocco  variables) 
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Equation  (A-32)  is  recognized  as  that  obtain.  ’ empirically  by  Danberg  (when  i is  unity). 


In  the  above,  it  is  noted  that  for  adialjatic  wall  conditions  (T  = T ),  is  zero 

aw  w 

and  one  obtains  the  parabolic  distribution  that  is  characteristic  of  nozzle  wall  data 
(for  Pr  unity).  On  the  other  hand,  when  considering  a unit  Prandtl  number  condi- 
tion where  T - T°  and  P is  unity,  Eq.  (A-32)  reduces  to 
aw  e j -1  X / 
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which  is  representative  of  the  classic  linear  Crocco  relation  (for  { = unity).  Thus  the 


important  contribution  of  Kutateladze  and  Leont’ev  is  the  parameter  I = ^ 


(f)' 
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Figure  1.  Model  Schematic 


Re^ 

Tw/C 

Model 

6.2 

1.7  X 10^-3.5  X 10^ 

.56- .84 

Flat  plati: 

6.4 

6x  10^ 

.52 

Flat  plate 

6.0 

1.1  X 10^ 

.38.  .49 

Hollow  cylinder 

5.1 

3.1  X 10^-4.0  X 10^ 

.72-83 

Flat  plate 

10.2 

2.3  X 10^ 

78 

Cone 

3.67 

9x10^ 

.42 

Parabolic,  finenett  ratio  of  10 

6.75 

3S  X 10^ 

.63 

Hollow  cylinder 

10.6 

1.3  X 10^ 

.3 

Flat  plate 

5,6.8 

2x  103-1.3x  10^ 

.4-.7 

Hollow  cylinder 

6.6 

2.2  X 10^-6.9  X 10^ 

1 

.3-38 

1 

Flat  plate 

0 I .2  .3  .4  .6  .6  7 8 9 10 


Figure  2.  Total  Temperature  - Volooi^  for  "Flat  Plate"  Flows  3.67  <.  10..'); 

0.28  T.JT^  <s  0.84  (Ref.  6) 
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Figure  9.  Profile  Data  Summary  for  a Transitional-Turbulent  Boundary  Layer 
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Figure  47.  Empirical  Heat  Blockage  Correlations  of  Previous  Investigators 
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Verification  of  the  Heat  Blockage  Correlation  (Foreign  Gas  Injection) 
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Figure:  r.4.  Comparison  of  Heat  Blockage  Correlation  with  Air  Injection  Data 
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